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ABSTRACT 

The detailed abundances of 23 chemical elements in nine bright RGB stars in the Carina dwarf 
spheroidal galaxy are presented based on high resolution spectra gathered at the VLT and Magellan 
telescopes. A spherical model atmospheres analysis is applied using standard methods (LTE and 
plane parallel radiative transfer) to spectra ranging from 380 to 680 nm. Stellar parameters are 
found to be consistent between photometric and spectroscopic analyses, both at moderate and high 
resolution. The stars in this analysis range in metallicity from — 2.9 < [Fe/H] < —1.3, and adopting 
the ages determined by Lemasle et al. (2012), we are able to examine the chemical evolution of Carina's 
old and intermediate-aged populations. One of the main results from this work is the evidence for 
inhomogeneous mixing in Carina; a large dispersion in [Mg/Fe] indicates poor mixing in the old 
population, an offset in the [a/Fe] ratios between the old and intermediate-aged populations (when 
examined with previously published results) suggests that the second star formation event occurred 
in a-enriched gas, and one star, Car-612, seems to have formed in a pocket enhanced in SN la/II 
products. This latter star provides the first direct link between the formation of stars with enhanced 
SN la/II ratios in dwarf galaxies to those found in the outer Galactic halo (Ivans et al. 2003). Another 
important result is the potential evidence for SN II driven winds. We show that the very metal-poor 
stars in Carina have not been enhanced in AGB or SN la products, and therefore their very low ratios 
of [Sr/Ba] suggests the loss of contributions from the early SNe II. Low ratios of [Na/Fe], [Mn/Fe], 
and [Cr/Fe] in two of these stars support this scenario, with additional evidence from the low [Zn/Fe] 
upper limit for one star. It is interesting that the chemistry of the metal-poor stars in Carina is not 
similar to those in the Galaxy, most of the other dSphs, or the UFDs, and suggests that Carina may 
be at the critical mass where some chemical enrichments are lost through SN II driven winds. 
Subject headings: galaxies: abundances — galaxies: evolution — galaxies: dwarf — galaxies: individ- 
ual (Carina dwarf galaxy) 



1. INTRODUCTION 

Chemical analyses of stars in nearby classical dwarf 
galaxies have swelled in the past decade, with the ad- 
vent of large aperture telescopes, high efficiency spec- 
trographs, and multiplexing capabilities. Detailed abun- 
dances are now available for dozens to hundreds of stars 
in classical dwarf spheroidal galaxies (dSph) , ultra faint 
dwarf (UFD) systems, and even a few massive stars in 
some low mass dwarf irregular galaxies (dirr); see the 
review by Tolstoy et al. (2009). These chemical stud- 
ies have shown that low mass dwarf galaxies have had 
a slower chemical evolution than the stellar populations 
in the Milky Way. For example, the majority of stars 

[email: kvenn@uvic.ca 

* This project is partially based on VLT FLAMES spectroscopic 
observations obtained at the European Southern Observatory, 
proposals 074.B-0415 and 076.B-0146, and partially based on 
MIKE spectra gathered at the 6.5-m Magellan Telescopes, in 
Chile. 



in dSphs have mean metallicities and [a/FeQ ratios that 
are lower than those of the Sun. Stars in dwarf galax- 
ies also tend to have different heavy element ratios (e.g., 
higher [Ba/Y] or [La/Eu]) than simila r metallicity stars 
in the Galaxy (e.g., Venn et al. 2004, lAoki et al. 1120091 
Letarte et al. 2010). On the other hand, the most metal- 
poor stars in dwarf galaxies, with [Fe/H] < —2.5, tend 
to have similar abundance ratios to Galactic halo stars, 
even for the a and heavy element ratios (e.g., Sculptor, 
Fornax, and Sextans, Tafelmeyer et al. 2010). This is 
also true of the very metal-poor stars in the UFDs (e.g.. 
Com Ber, Bootes I, and Leo IV, Norris et al. 2008, Feltz- 
ing et al.2009, Frebel et al. 2010a, Simon et al.2010). 
One exception are the a and heavy element abundances 
for metal-poor stars in Sextans, which show both offsets 
and larg er dispersions tha n the Galaxy and other dwarf 
galaxies (|Aoki et al. 1120091 ). 

1 [X/Fe] = log(X/Fe). - log(X/Fe)0. 
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The Carina dwarf galaxy provides a new opportunity 
for chemical evolution and nucleosynthetic studies. It 
has a low mass that is similar to that of Sextans (Walker 
et al. 2009), and it has had an unusual, episodic star for- 
mation history. Its colour-magnitude diagram (CMD) 
shows at least three main sequence turn-offs (Monelli 
et al. 2003, Bono et al. 2010) and is best described by a 
star formation history with a well-defined old population 
(^^10-12 Gyr), a dominant intermediate-aged population 
(~5-7 Gyr), and a trace young population (~2 Gyr; al- 
though the specific timescales are uncertain). These star 
forming episodes are separated by long quiescent periods, 
seen as gaps between the main sequence turn-off points. 
It is estimated that 70 to 80% of the stars in Carina have 
intermediate- ages, while most of the remaining stars are 
old and associated with the first star forming episode 
(also see Dolphin 2002, Hernandez et al. 2000, Hurley- 
Keller et al. 1998, and Migheh 1997). One wonders how 
this unique star formation history may have affected the 
chemical evolution of Carina. 

In spite of its punctuated star formation history. Ca- 
rina has an extremely narrow red giant branch. One pos- 
sibly is that the RGB stars in Carina have a fortuitous 
alignment in the age-metallicity degeneracy, such that 
the older metal-poor stars overlay the metal enhanced 
intermediate-aged stars. Alternatively, the narrow red 
giant branch could be dominated by intermediate- aged 
stars only, if that population formed on a relatively short 
timescale and with only a modest metallicity spread 
(Rizzi et al. 2003, Bono et al. 2010). Low resolution spec- 
tra of 437 red giant branch (RGB) stars in Carina anal- 
ysed by Koch et al. (2006) showed the mean metallicity 
in Carina is [Fe/H] ~ — 1.7± 0.9, in agreement with sim- 
ilar analyses by Helmi et al. (2006) and Starkenburg et 
al. (2010). This meta l hcity spread is larger than pre- 
dicted by iBono et al.l (|2010( ) from their CMD analysis 
(A[Fe/H] < ±0.5 dex). Koch et al. (2006) also found no 
significant gradients in metallicity or stellar population 
with position in Carina, unlike the results for Sculptor 
(Tolstoy et al. 2004) and Fornax (Battaglia et al.2006). 

High resolution analyses of five RGB stars (Shetrone 
et al. 2003) and ten RGB stars (Koch et al. 2008) sup- 
ported the larger range in metallicity found by the low 
resolution spectral analyses. These analyses also showed 
that Carina is like the other dSph galaxies in that the 
[a/Fe] ratios are lower than in Galactic halo stars with 
similar metahicity. Lanfranchi et al. (2003, 2004, 2006) 
used these datasets to develop chemical evolution mod- 
els for Carina, tuning the high wind efficiency and low 
star formation efficiency to reproduce Carina's metallic- 
ity distribution function, [a/Fe] ratios, and low gas con- 
tent. 

In this paper, we present a new detailed abundance 
analysis of up to 23 elements in nine stars in the Carina 
dwarf galaxy. This work increases the number of ele- 
ments, the number of stars, and the metallicity range 
previously explored. We also report the abundances 
for two newly discovered very metal-poor stars (with 
[Fe/H]^ —2.85), which allow us to examine the earliest 
epoch of star formation. These results, with previously 
published abundances, are used to explore the unique 
star formation history and chemical evolution in Carina. 

2. OBSERVATIONS AND DATA REDUCTIONS 
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Fig. 1. — Colour-magnitude diagram, V vs (V-I), for Carina from 
our ESO WFI data. Only stars with radial velocities > 200 kms~^ 
(from Koch et al. 2006) are shown. Stars with high resolution 
spectral analyses are identified by red symbols, including nine 
stars from this paper, five from Shetrone et al. (2003), and ten 
from Koch et al. (2008; four are in common with this analysis). 
One carbon-star is noted by the red cross, and two foreground 
objects are shown as black crosses. The faintest star in our sample 
is Car-7002. 



The data presented in this paper were acquired at two 
observatories during separate time allocations in January 
2005. The FLAMES (Fibre Large Array Muhi Element 
Spectrograph) multiobject spectrograph at the 8.2 me- 
ter UT2 (Kucyen) at the Very Large Telescope (VLT) 
of the European Southern Observatory (ESO) was used 
to collect high resolution spectra using both the UVES 
and GIRAFFE fiber modes (Pasquini et al.2002). The 
analysis of the FLAMES/GIRAFFE spectra is presented 
by Lemasle et al. (2012). Also, the the MIKE (Magellan 
Inamori Kyocera Echelle) spectrograph at the Magellan 
Landon 6.5m Clay Telescope at the Las Campanas Ob- 
servatory was used to collect high resolution spectra of 
individual stars outside of the central field of Carina. Our 
targets and those previously analysed with high resolu- 
tion spectral analyses by Koch et al. (2008) and Shetrone 
et al. (2003) are shown in Fig.[T]on the V vs (V-I) CMD. 
These stars are also shown in Fig. [2] as metallicity (from 
Koch et al. 2006) versus their location within Carina. 

2.1. Magellan MIKE Spectra 

Target selection for the Magellan run was from low 
resolution spectra of the Ca II 8600 A feature (=CaT) 
for hundreds of stars in the Carina dSph in the ESO 
archive (Koch et al. 2006, Helmi et al. 2006). The feature 
was used to ascertain membership from radial velocities, 
as well as estimate an initial metallicity for each tar- 
get based o n the CaT-metallicity calibration available at 
that time (jBattaglia et al.l 120081 ). Exposure times, ra- 
dial velocities (kms~^), the signal-to-noise ratio (SNR) 
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Fig. 2. — Spatial location and metallicities of the nine Carina 
stars in this paper, as well as five from Shetrone et al. (2003), and 
ten stars from Koch et al. (2008; four are in common with this anal- 
ysis). Symbols are the same as in Figure [l] These are compared 
to the distribution of 437 RGB stars from Koch et al. (2006) which 
have a high (3 cr) probability of Carina membership. The metallic- 
ities on this plot are those determined by Koch et al. (2006). Ellip- 
tical radii are calculated using the structural parameters for Carina 
from Irwin & Hatzidimitrou (1995; central coordinates 6''40.6' and 
-50°56' for epoch B1950, position angle 65°, ellipticity e=0.33), 
and listed in Table [T] and core and tidal radii are indicated. 



at three wavelengths, the elliptical radius from the cen- 
ter of Carina (in arcminutes), and alternative names for 
the targets are summarized in Table [TJ Spectra for three 
objects in the metal-poor globular cluster M68 were also 
taken as standard stars. M68 was chosen as a standard 
since it has low metallicity (Harris 1996) and low red- 
dening. BVI magnitudes are from Stetson (2000, with 
corrections onto the Johnson Kron- Cousins scale of the 
Landolt standard stars as in Stetson 2005), which places 
these stars on the upper RGB of M68 and suggests they 
have intrinsic luminosities similar to our program stars. 

On the red side, the double echelle design covers 4850 - 
9400 A, and on the blue side, 3800 -5050 A. The quality of 
the blue spectra is significantly higher in the overlapping 
wavelength region (~ 4850 to 5050 A). Using a 1.0" slit, a 
spectral resolution R = 28,000 (blue) to R = 22,000 (red) 
was obtained. The red chip had a gain of 1.0 electrons 
ADU~^ and read noise of 4 electrons; the blue chip had 
a gain of 0.47 electrons ADU~^ and a read noise of 2 
electrons. We binned on-chip with 2x2 pixels. Table [T] 
lists the SNR per pixel achieved in the final combined 
spectra near 4200, 5200, and 6200 A. 

The data were reduced using standard IRAjQ routines. 
Sky subtraction was done with a smooth fit perpendicu- 
lar to the dispersion axis. Heliocentric corrections were 
applied to each spectrum before determining their ra- 
dial velocities. Multiple spectra were median combined, 
which helped to remove cosmic ray strikes. Spectra taken 
of the Th-Ar lamp provided the wavelength calibration. 
Both quartz flats and screen flats were taken at vary- 
ing exposure levels; no significant offsets were found be- 
tween the well-illuminated quartz flats and the science 
exposures, thus the quartz flats were adopted for the 
data reductions. A small dark current was noticed on 
one corner of the blue chip in a region that did not re- 
ceive starlight and therefore did not affect this analysis. 
The final spectra were normalized using k-sigma clipping 
with a non-linear filter (a combination of a median and 

^ IRAF is distributed by the National Optical Astronomy Ob- 
servatory, which is operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 



a boxcar). The effective scale length of the filter was set 
from 8 to 15 A, dependent on the crowding of the spec- 
tral lines, and we found that this was sufficient to follow 
the continuum without affecting the presence of the lines 
when used in coniunctio n. with iterative c lippi ng. This 
method was also used by iBattaglia et al. I ()2008l) to nor- 
malize CaT spectra. 



2.2. VLT FLAMES/UVES Data 

Observations with the multi-object FLAMES spec- 
trograph at the UT2 Kueyen ESO-VLT (Pasquini et 
al. 2002) were carried out for nine targets in the cen- 
tral region of Carina. In two separate configurations, 
the FLAMES/UVES fibers were placed on bright RGB 
stars resulting in nine stellar spectra and two sky ob- 
servations. Two of the nine targets proved to be fore- 
ground RGB stars, while one is a carbon star that is un- 
suitable for our analysis. Target coordinates, exposure 
times, radial velocities (kms~^), elliptical radius from 
the center of Carina (in arcmins), and the SNR at two 
wavelengths for the final combined spectra are listed in 
Table [TJ Other names for each target are also listed. 
FLAMES/GIRAFFE spectra were simultaneously ob- 
tained at high resolution (R=20,000) for 36 more stars 
in Carina, but over much shorter wavelength regions (3 
wavelength settings that yielded ~250 A each) - these are 
presented in a separate paper by Lemasle et al. (2012). 
Due to variable weather conditions, only eleven of the 
one-hour exposures had sufficient signal to be used fur- 
ther in this analysis. 

On the red side, the double echelle design covers 
5840-6815 A, and on the blue, 4800-5760A. A 1.0" slit 
yields a spectral resolution R = 47,000. Table [T] lists the 
exposure times and SNR per pixel in the final combined 
spectra near 5200 and 6200 A. Unfortunately, six of the 
one-hour exposures did not have enough signal in the 
FLAMES/UVES fibers to be useful in this analysis. The 
data were reduced using the ESO FLAMES pipeline0 
The sky spectrum was fit with a smoothly varying func- 
tion, and this was subtracted from the stellar exposures 
per wavelength set up (to reduce adding more noise in 
already low SNR spectra). Each spectrum was then he- 
liocentric velocity corrected, radial velocities were deter- 
mined, and the spectra were coadded (weighted by the 
SNR). The final spectra were normalized using an iter- 
ative asymmetric k-sigma clipping routine with a non- 
linear filter (like for the Magellan/MIKE spectra). 

3. STELLAR ANALYSES 

3.1. Line List 

A range of elements are detectable in our spec- 
tra, which enables a comprehensive abundance analy- 
sis. Atomic lines for this analysis were selected from 
the literature, including the lin e lists from Shetrone et 
al. (2003), Cayrel et al. (2004), lAoki et al.l (l2007l) . Co- 
hen et al. (2008), Letarte et al. (2009), Tafelmeyer et 

3 The GIRAFFE Base-Line Data Reduction, girBLDRS, 
was written at the Observatory of Geneva by A. Blecha 
and G . Simond. and is ava ilable through SourceForge at 
|http: / /gir bldrs. sourceforge.net/l 
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TABLE 1 

Observing Information 
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Note. — 


Data for M68 and the first three Carina stars 


are from the Magellan MIKE spectrograph, with the remaining 


nine 



Carina stars observed with the FLAMES spectrograph at the VLT. Stars partially analysed by Koch et al. (2008) are noted 
with "(K)" beside their alternative names. The S# target names are from P.B. Stetson's onhne database of homogeneous 
photometry (at http://cadcwww.hia.nrc.ca/stetson). Car-X is a bright target right next to Car-909, which was the actual 
target. SNR are the signaFto-noise^raSos per p%xeL 

al. (2010), and Frebel et al. (2010a); see TableEl Atomic 
data for the Fe I lines were updated from O 'Brian et 
al. (1991) when available, or the atomic data was up- 
dated to the latest values in the National Institute of 
Standards and Technology (NIST) databaseQ 



3.2. Equivalent Width Measurements 

Most of the elemental abundances in this analysis are 
determined from equivalent width (EW) measurements. 
Spectrum syntheses were used only for lines affected by 
hyperfine splitting or for elements with line measure- 
ments from low SNR spectra. 

Equivalent widths were measured using the Gaussian 
fitting routine DAOSPEC (Stetson & Pancino 2008). 
The placement of the continuum is a critical step, as 
it influences the measurements of the equivalent widths; 
thus some care was taken in assigning and testing the 
DAOSPEC fitting parameters. For example, a low order 
polynomial (order 13) was adopted to allow DAOSPEC 
to measure the effective continuum, and tests of other 
low order values (order =5 and =1) resulted in nearly 
identical equivalent widths, with AEW < 1 %. As a final 
exercise, the continuum was set to one so that DAOSPEC 
would not relocate the global continuum for the measure- 
ment (order = -1). Results from this test showed a slight 
offset < -|-12% in the equivalent width measurements, as 
expected. We note this offset is similar to or less than our 

http:/ /physics. nist.gov/PhysRefData/ASD/index. html 



adopted equivalent width measurement errors from the 
Cayrel formula (described below). For testing purposes, 
equivalent widths were also measured with splot in IRAF 
to determine the area under the continuum. A compar- 
ison of these measurements with those from DAOSPEC 
for two of the stars observed with the MIKE spectro- 
graph are shown in Fig. [3l 

The expected minimum measurement uncertainties on 
the equivalent widths (EWrms) have been estimated us- 
ing a revised Cayrel formula, i.e., a new derivation of the 
measurement error s on th e equiv alent widths was pre- 
sented by Battagli a et al.l (|2008[ ). where it was found 
that the factor of 1.5 in the Cayrel (1988) formula is 
actually within the square root; thus, 

EWrms = SNR-i x Vl.5 X FWHM x Sx 

where SNR is the signal-to-noise ratio per pixel, FWHM 
is the line full width at half maximum, and 5x is the pixel 
size. An extra 10% of the EW was added to this for a 
more conservative measurement error, such that 

AEW = EWrms + 0.10 X EW. 

For the M68 stellar spectra taken at Magellan the re- 
vised Cayrel formula gives EWrms = 3 mA, but EWrms — 
6 mA for the lower SNR spectra of the Carina stars ob- 
served with MIKE. For the spectra taken at the VLT 
with FLAMES/UVES, the higher resolution yields a 
minimum error of EWrms = 4 mA, with an exception 
for Car-1013 where EWrms = 7 mA due to its lower SNR 
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TABLE 2 
Line List 



Wave X log gf CAR CAR CAR CAR CAR CAR 
(A) (eV) 484 524 612 705 769 1013 



Fe I 

5001.864 3.88 0.010 124.9 116.7 124.0 134.0 104.3 165.1 

5006.120 2.83 -0.615 194.7 166.1 185.2 193.8 158.2 

5012.070 0.86 -2.642 •■• 195.0 179.3 

5014.940 3.94 -0.300 113.4 103.4 117.0 114.3 65.4 97.2 

Note. — Equivalent widths are in mA, and upper limits are noted. 
When spectrum syntheses have been used in the abundance analysis, then 
an "S" is noted. A second table includes the stars observed with the MIKE 
spectrograph at Magellan. Both tables are published in their entirety in 
the electronic edition. A portion is shown here for guidance regarding its 
form and content. 



spectrum. 

The reported errors from the DAOSPEC program were 
lower than those from the revised Cayrel formula, by 
1/3 to 1/2. This is due to correlations in the noise 
estimates in DAOSPEC when the pixel data are re- 
binned/interpolated during the spectral extractions and 
wavelength calibrations (P. Stetson, priv. communica- 
tions). When individual line measurements were signif- 
icantly different between the DAOSPEC and splot mea- 
surements, those lines were checked by eye. These lines 
tended to be unresolved blends and / or suffered from dif- 
ficult continuum placement. Since DAOSPEC uses a 
fixed FWHM and a consistent prescription for contin- 
uum placement, those EW measurements were adopted 
for this analysis, although the AEW errors are from the 
revised Cayrel formula. 




50 100 150 200 

DAOSPEC EW 



Fig. 3. — Equivalent width measurement comparisons for M68- 
6022 and Car-1087, both observed with the MIKE spectrograph. 
The three lines represent equal EW values and EW it AEW, where 
AEW = EWrms -I- 10% EW, and EWrms = 3 mA and 6 mA, 
respectively. 



A maximum EW of 200 mA was adopted for all 
stars since deviations from the Gaussian profiles used 
in DAOSPEC become more significant for stronger lines. 
In fact, for lines >200 mA, we found that DAOSPEC 
would occassionally divide these lines into 2-3 lines, each 
with lower EWs. The final line lists and measurements 
were carefully reviewed during this analysis. 

3.3. Radial Velocities 

Radial velocities were initially measured from 5-7 
strong lines and refined using DAOSPEC (which cross 
correlates the measured positions of the detected spectral 
lines with reference wavelengths from an input linelist). 
Heliocentric corrections were then applied. The radial 
velocities for the Carina stars are listed in Table [1] the 
results for the Magellan/MIKE blue and red spectra were 
averaged together since each arm was reduced indepen- 
dently. The three M68 stars have a mean radial velocity 
of — 95.4± 1.0 kms~^, with an average measurement er- 
ror of 1.5 kms~^. This is in agreement with results in 
the literature, e.g.. Lane et al. (2009) report a radial ve- 
locity of -94.93 ± 0.26 kms'^ from 123 RGB stars in 
M68 taken at the 3.9m Anglo Australian Telescope with 
the AAOmega spectrograph. The mean radial velocity 
of our nine Carina stars is 220.5 ± 4.9 kms~^; the aver- 
age measurement error of the FLAMES /UVES data are 
0.9 kms~^, whereas it is 3.8 kms"^ for the metal-poor 
stars with lower quality Magellan/MIKE spectra. These 
results are in agreement with the mean radial velocity of 
223.9 kms-i found by Koch et al. (2006) for 437 RGB 
members of Carina, with a dispersion of 7.5 kms~^. The 
radial velocities from both the high resolution (this pa- 
per) and lower resolution spectra (from Koch et al. ) are 
shown in Table [TJ again, these are in good agreement, 
with a mean difference in ( RVcaT — IWhr ) = 2.1 ± 
3.4 kms"^ 

3.4. Photometric Parameters (physical gravity) 

BVIJKs photometric values for the Carina targets and 
M68 standard stars are listed in Table \3\ BVI magni- 
tudes for the M68 stars are from Stetson (2000, on the 
Johnson Kron-Cousins scale. Stetson 2005). J and Kg 
magnitudes are from the 2MASS database0. For the Ca- 
rina members, VI magnitudes are from the ESO 2.2m 

^ This publication makes use of data products from the Two 
Micron All Sky Survey which is a joint project of the Uni- 
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TABLE 3 
Target Information 



Star BVIVI J KsBVI JK 

PBS PBS PBS WFI WFI 2MASS 2MASS Gul Gul Gul Gul Gul 



M68-6022 


15.127 


14.205 


13.106 






12.324 


11.659 












M68-6023 


15.192 


14.284 


13.191 






12.363 


11.712 












M68-6024 


15.230 


14.354 


13.289 






12.479 


11.858 












Car-1087 


19.146 


18.031 


16.735 


18.00 


16.73 


15.802 


15.301 












Car-5070 


19.125 


17.904 


16.520 


17.86 


16.49 


15.641 


14.758 












Car-7002 


19.284 


18.344 


17.146 


18.40 


17.19 


16.244 


15.562 












Car-484 


19.009 


17.603 


16.179 


17.56 


16.16 


15.322 


14.627 


19.055 


17.607 


16.178 


15.290 


14.406 


Car-524 


18.934 


17.645 


16.308 


17.63 


16.32 


15.388 


14.473 


18.921 


17.598 


16.278 


15.398 


14.574 


Car-612 


19.077 


17.811 


16.416 


17.82 


16.42 


15.456 


14.763 


19.103 


17.779 


16.403 


15.609 


14.828 


Car-705 


19.101 


17.828 


16.430 


17.78 


16.43 


15.687 


14.597 


19.104 


17.788 


16.424 


15.509 


14.663 


Car-769 


19.202 


18.002 


16.693 


18.01 


16.71 


15.698 


14.921 


19.224 


18.025 


16.725 


15.758 


15.040 


Car-1013 


19.048 


17.783 


16.424 


17.77 


16.42 


15.515 


14.888 


19.056 


17.767 


16.421 


15.528 


14.726 


Car-489 








17.85 


16.53 


15.544 


14.933 












Car-837 








17.71 


16.47 


15.456 


14.509 


19.111 


17.705 


16.509 


15.427 


14.621 


Car-X 








17.51 


16.42 


15.748 


15.223 













Note. — BVI from P.B. Stetson's website of homogeneous photometry, VI from our WFI data, JK^ from the 2MASS 
database, and BVIJK from M. GuUeuszik (private communications). 



WFI photometry; initial ealibrations were on the defauh 
ESO zero-point values, but have been updated to the 
Johnson Kron-Cousins scale with Stetson's database and 
Gullieuszik's photometry (private communications). 

Initial metallicities for the Carina stars are from the 
CaT measurements (Koch et al. 2006). For the three M68 
stars, the average metallicity from Lee et al. (2005) is 
adopted, [Fe/H]=-2.16. 0. 

Color temperatures were found using the Ramirez & 
Melendez (2005) calibration, and adopted as the initial 
effective temperatures; see Tabled Deviations between 
the different color temperatures are quite small for the 
three M68 globular cluster stars, i.e., ctt < 20 K from 
T{B - V), T{V - I), T{V - J), and T{V - K). Devia- 
tions are larger for the Carina stars, with ctt < 134 K, 
and ((Tt) — 74 K. Physical (photometric) gravities are 
determined using the standard relation. 



log g = log Qq 



log (nt) + 4xlo; 
+ 0.4x(Mbo1: 



MboIo) 



with log 50 = 4.44, Teff© = 5777 K, MboIq = 4.75, and 
adopting a reddening law of A(V)/E(B-V) = 3.24. The 
values of MboI* are based on the bolom etric correction 
to the V magnitude from lAlonso et al.l ([l999^ , with an 
assumed distance and stellar mass. 



reddening is adopted. A turn off mass of 0.83 ±0.03 Mq 
is found for M68 using isochrones by Vandenberg & Bell 
(1985). 

For Carina, Mighell (1997) determined the distance 
modulus and reddening from WFPC2 V and I band 
imaging, but also compares his results in an appendix 
to the many estimates available in the literature. We 
adopt Mighell's determination of (m-M)v=20.05 ±0.11 
and E(B-V)=0.06 ±0.02. This reddening value is the 
same as from the Schlegel et al. (1998) maps. One of 
the highest values for the distance modulus was found 
by Monelh et al. (2003; (m-M)v = 20.24, with E(B- 
V)=0.03), but as discussed below (also see Section 5); 
these differences have negligible effects on our spectral 
analysis. An estimate of the turn off mass is more diffi- 
cult in Carina due to its separate episodes of star forma- 
tion, each of which has a different turn off mass; these 
estimates can have a larger effect on the physical gravity 
than the distance modulus. Since Carina is dominated 
by an intermediate aged population (5-7 Gyr old), we ex- 
amined isochrones (Fagotto et al. 1994, with Y=0.23 and 
Z=0.0004) to find that the age range of 5 to 12 Gyr corre- 
sponds to turn off masses of 1.0 to 0.8 M©, respectively. 
We adopted 0.8 M© in this analysis, but note that if 1.0 
Mq were used, this would change the physical gravities 
by +0.10 dex with no effect on temperature; this change 
has only a small effect on the chemical abundances. 



For M68, we adopt the distance and reddening in the 
Harris catalogue (1996), (m-M)v = 15.19 ±0.10 and E(B- 
V)=0.05 from isochrone fitting by McClu re et al. (1987) 
This value is in good agreement with IBrocato et al.l 
(|1997|1 . but the proper motion distance modulus, (m- 
M)v = 14.91 (Dinescu et al. 1999) is quite short if this 

versify of Massachusetts and the Infrared Processing and Anal- 
ysis Center/California Institute of Technology, funded by the 
National Aeronautics and Space Administration and the Na- 
tional Science Foundation. The database can be found at 
http: / / www.ipac.caltech.edu/2mass/releases / allsky 

" Lee et al. (2005) found a mean [Fe II/HJ=— 2.16 from a high 
resolution spectral analysis of eight stars in M68, using photomet- 
ric gravities. They also find a mean [Fe l/H]=— 2.56 using spectro- 
scopic gravities. 



3.5. Spectroscopic Parameters 

The CaT metallicity and the photometric parameters 
for temperature and gravity were adopted as the ini- 
tial stellar parameters for a model atmospheres analysis. 
Spectroscopic indicators were used to refine the effective 
temperature and metallicity, as well as to determine the 
microturbulcnce values. These values are listed in Ta- 
ble [5l Only the physical gravity is unchanged from the 
photometric parameters analysis. 

3.5.1. Model Atmospheres 

The new MARCS spherical models have been adopted 
in this analysis (Gustafsson et al. 2003, 2008, further ex- 
panded by B. Plez, private communications). The grid 
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TABLE 4 

Photometric Parameters 



Target 



Tbv 
(K) 



Tvi 

(K) 



Tvj 
(K) 



Tvk 
(K) 



(T) 
(K) 



(K) 



Mv Mbol 



M68-6022 

M68-6023 

M68-6024 

Car-1087 

Car-5070 

Car-7002 

Car-484 

Car-524 

Car-612 

Car-705 

Car-769 

Car-1013 



4663 
4686 
4741 
4391 
4303 
4677 
4073 
4228 
4229 
4194 
4308 
4208 



4671 
4682 
4732 
4457 
4314 
4552 
4245 
4346 
4253 
4305 
4361 
4280 



4677 
4654 
4722 
4572 
4367 
4492 
4383 
4235 
4301 
4221 
4278 
4372 



4694 
4647 
4701 
4416 
4364 
4473 
4320 
4311 
4224 
4460 
4254 
4250 



4676 
4667 
4724 
4459 
4337 
4548 
4255 
4280 
4252 
4295 
4300 
4278 



13 
20 
17 
80 
33 
92 

134 
58 
35 

120 
46 
70 



1.53 
1.55 
1.61 
0.95 
0.83 
1.16 
0.65 
0.70 
0.75 
0.77 
0.86 
0.75 



-0.98 
-0.91 
-0.84 
-2.05 
-2.19 
-1.65 
-2.49 
-2.42 
-2.23 
-2.27 
-2.04 
-2.28 



-1.40 
-1.33 
-1.23 
-2.59 
-2.79 
-2.15 
-3.14 
-3.05 
-2.88 
-2.89 
-2.65 
-2.91 



Note. — Photometric temperatures are from the Ramirez & Melen- 
dez (2005) calibration. Each temperature is determined from a photometric 
colour, e.g., Tbv represents the temperature determined from the (B— V) 
colour. The average of the four colour temperatures ((T)) and standard 
deviation (ctt) are hsted along with other parameters derived from the mag- 
nitudes (see Section 13.41 1. 



covers the range of parameters (temperature, luminos- 
ity, microturbulence, and metallicity) of RGB stars. The 
models used in this analysis adopt the Galactic abun- 
dance pattern, i.e., [a/Fe] = +0A for [Fe/H] < -1.0, 
which is the metallicity range of the stars in Carina. We 
adopt the standard version of the MOOGd (Sneden 1973) 
spectral synthesis code to determine the individual line 
abundances. Corrections were applied when necessary 
for continuum scattering effects that can affect the re- 
sults for lines below 5000 A, as described in Section [3.8.2l 
The photometric gravities were retained throughout the 
analysis (see Section [3.8.3p . 



3.5.2. Effective Temperatures, Microturbulence, 
Metallicity 



& 



The effective temperatures were refined from the pho- 
tometric values through examination of the Fe I lines. 
Only lines redwards 5000 A were considered to avoid un- 
certainties in the standard MOOG treatment of scatter- 
ing in the continuous opacity (see Section I3.8.2P . Micro- 
turbulence values were found by minimizing the slope in 
the Fe I line abundances with equivalent width. Simi- 
larly, the excitation temperatures were found by mini- 
mizing the slope in the same Fe I line abundances with 
excitation potential {x, in eV). While the M68 stars had 
lower excitation than colour temperatures (by an average 
of 67 K), the Carina stars had higher excitation tempera- 
tures (by an average of 200 K). The final metallicity (for 
which we adopt the [Fe I/H] values) is also found at the 
end of these minimization iterations. 

No significant relationship was found in the Fe I abun- 
dances versus wavelength (>5000 A) for either the glob- 
ular cluster or the Carina targets; this indicates that the 
sky subtraction and overall data reduction methods were 
successful. Eliminating lines with x < 1.4 eV, which can 
be particularly sensitive to departures from local thermo- 
dynamic equilibrium (LTE) effects, also had no signifi- 
cant effect on the temperature or microturbulence values. 



^ T he 2009 version of MOOG i s 
|http: / /www.as.utexas.edu/'-'.^chris/moog.htmH 



available from 



Samples of these relationships are shown in Fig.|4]for two 
stars, one observed with the FLAMES/UVES spectro- 
graph and the other observed with the MIKE spectro- 
graph. Adjusting the metallicities from the initial CaT 
values had no significant effect on the excitation tempera- 
tures (AT <10 K). Adopting the spectroscopic tempera- 
tures had only a small effect on the photometric gravities 
(< +0.08 in the Carina stars, and < -0.02 for the M68 
stars). 

Four stars in Carina have been observed and analysed 
independently by Koch et al. (2008). Their analysis with 
ATLAS9 models atmospheres results in excitation tem- 
peratures that are in excellent agreement with ours, with 
differences ranging from A T = -1-116 to -70 K, and with 
an average offset of -1-32 K. The differences in microtur- 
bulence are more significant, with Koch et al. 's values 
being ~ 0.50 kms~^ lower. This is linked to differences 
in our gravity determinations discussed in Section [3.8.31 

3.6. Stellar Parameter Error Estimates 

Stellar parameters and their uncertainties are listed 
in Table [5l The one sigma uncertainty in the effective 
temperature is determined from the slope in the Fe I 
line abundances vs excitation potential (%) allowing the 
change in the iron abundance to equal the standard de- 
viation a(Fe I. Similarly, the one sigma uncertainty in 
microturbulence is determined by setting the slope of the 
Fe I line abundances vs equivalent width equal to tT(Fe I). 

For the physical gravities, the uncertainties are deter- 
mined from errors in the distance moduli, reddening, and 
stellar mass (see Section (23) . For M68, the uncertainty 
in the turn off mass is tiny (±0.03 Mq), which has a neg- 
ligible effect on the gravity (Alog5= ±0.04). The uncer- 
tainties in gravity are also very small when the short 
distance modulus reported by Dinescu et al. (1999) is 
adopted for M68 (Alog5= ±0.07). For Carina, the dom- 
inant uncertainty is from the stellar mass. The range in 
turn off masses for Carina is 0.8 to 1.0 Mq due to its 
episodic star formation history; the resulting uncertain- 
ties in log 5 can be as large as 0.10 dex. 

Uncertainties in metallicity are adopted from the error 
in the mean of the Fe I abundances (see Section 14. 2p . 
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TABLE 5 

Spectroscopic Parameters (photometric gravities) 



Target 


Tcft- 


log 9* 




12+log(Fc I) 


[Fe I/H] ±cr (#) CaT 




(K) 


(cgs) 


(km/s) 


±5(Fo)=' 


(SIO)*^ 



M68-6022 

M68-6023 

M68-6024 

Car-1087 

Car-5070 

Car-7002 

Car-484 

Car-524 

Car-612 

Car-705 

Car-769 



4550 ±44 
4667 ±55 
4650 ±56 
4600 ±98 
4550 ±81 
4548 ±96 
4400 ±66 
4500 ±64 
4500 ±61 
4500 ±61 
4600 ±76 



1.53 ±0.04 
1.55 ±0.04 
1.61 ±0.04 
0.95 ±0.10 
0.83 ±0.10 
1.16 ±0.10 
0.65 ±0.10 
0.70 ±0.10 
0.75 ±0.10 
0.77 ±0.10 
0.86 ±0.10 



1.50 ±0.12 
1.60 ±0.12 
1.80 ±0.12 
2.45 ±0.25 
2.30 ±0.20 
2.00 ±0.24 
2.30 ±0.07 
2.40 ±0.07 
2.10 ±0.08 
2.10 ±0.07 
2.30 ±0.08 



4.99 ±0.02 
5.09 ±0.03 
4.86 ±0.03 
4.69 ±0.06 
5.35 ±0.05 
4.64 ±0.05 
5.97 ±0.02 
5.75 ±0.02 
6.20 ±0.02 
6.15 ±0.02 
5.82 ±0.02 



Car-1013 4600 ±105 0.75 ±0.10 2.20 ±0.11 6.20 ±0.03 



2.51 ±0.23 (92) 
2.41 ±0.24 (92) 
2.64 ±0.24 (81) 
2.81 ±0.34 (38) 
2.15 ±0.38 (55) 
2.86 ±0.33 (38) 
1.53 ±0.21 (114) 
1.75 ±0.21 (121) 
1.30 ±0.23 (126) 
1.35 ±0.25 (127) 
1.68 ±0.25 (125) 
1.30 ±0.37 (125) 



-2.16 
-2.16 
-2.16 
-3.10 
-2.57 
-3.29 
-1.62 
-1.61 
-1.84 
-1.63 
-1.63 
-1.58 



(5(Fe) is the error in the mean of the Fe I Unes (see Section 14.211 . 
^ CaT metaUicities from the caUbration by Starkenburg et al. (2010). For the three M68 stars, 
CaT metaUicities are from Lee et al. (2005). 
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Fig. 4. — MOOG results for Car-484 (left panel), a moderate SNR high resolution FLAMES/UVES spectrum, and Car-7002, a moderate 
SNR and moderate resolution Magellan/MIKE spectrum. These plots show that the Fe I line abundances are minimized with excitation 
potential (to determine surface temperature), equivalent width (to define microturbulence), and wavelength (a check on the sky subtraction). 



3.7. Spectrum Syntheses and Hyperfine Structure 
Corrections 

Spectrum syntheses were carried out for Na I (5688, 
5889, 5895), Mg I (3829, 3832), Sr II (4077, 4205), Ba II 
(4554, 5853, 6141, 6496), and Eu II (4129, 4205, 6645). 
In each case, the instrumental resolution of the spec- 
trum was adopted as the main source of broadening, and 
checked with the shapes and EW abundance results from 
nearby spectral lines (Fe, Ca, Ti). When an EW could be 
measured, the synthetic and EW abundance results were 
compared; when excellent agreement was found then the 
EW result was taken. Line abundances taken from the 
synthetic results are noted in Tabled 

Hyperfine structure (HFS) corrections were deter- 
mined from spectrum syntheses for the odd-Z elements 
(Sc II, V I, Mn I, Co I), Cu I, and the neutron cap- 
ture elements Ba II, La II, and Eu II. HFS components 
were taken from Lawler et al. (2001a; La II) , Booth 



et al. (1983; Mn I), Prochaska et al. (2001; Sc II, V I, 
Mn I, and Co I), and the Kurucz database0 for the re- 
maining lines. HFS c o mpon ents and isotopic ratios for 
Cu I are from iBiehfl ()1976f ). for Eu II from Lawler et 
al. (2001b), and for Ba II from McWiUiam et al. (1998; r- 
process isotopic ratios were adopted for the most metal- 
poor stars, whereas the solar ratios were used for stars 
with [Fe/H]> -2). 

3.8. Other Considerations 
3.8.1. Spherical vs Plane Parallel 

Spherical MARCS models represent a significant im- 
provement in the modeling of stellar atmospheres, e.g., 
Heiter & Eriksson (2006) recommend the use of spherical 
model atmospheres in abundance analyses for stars with 
\ogg< 2 and 4000 < Tcff< 6500 K, which encompasses 

« http://kurucz.harvard.edu/LINELISTS/GFHYPERALL 
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the range in stellar parameters of our target stars. Tests 
performed by Tafelmeyer et al. (2010) showed systematic 
offsets between the spherical and plane-parallel models 
are below 0.15 dex (this includes iron lines and other 
elements, e.g., [C/Fe]). 

We performed similar tests, comparing the metallic- 
ities we determined from spherical MARCS models to 
those from plane parallel MARCS (pp-MARCS) and Ku- 
rucz models. The plane parallel MARCS models had the 
effect of reducing the excitation temperature results by 
^100 K, but the Kurucz models had a different effect, 
to increase the microturbulence values by ^0.1 dex. In 
both cases, all other parameters were unchanged. These 
offsets, when applied to the abundance analysis, changed 
the Fe I abundances by only ^-^0.05 dex (+0.05 with pp- 
MARCS models, and -0.05 with Kurucz models). This 
is very similar to the offsets found by Tafelmeyer et 
al. (2010) and Heiter & Eriksson (2006) for iron. 

3.8.2. Continuum Scattering Effects 

The standard version of MOOG treats continuum scat- 
tering {ay) as if it were absorption {ny) in the source 
function, i.e., Sy = (the Planck function), an ap- 
proximation that is only valid at long wavelengths in 
RGB stars. At shorter wavelengths (<5000 A), Cayrel et 
al. (2004) have shown that the scattering term must be 
taken into account such that the source function becomes 
Sy = {ny By + ay J y) / + d y) . A ucw vcrslou of MOOG 
(MOOG-SCAT, Sobeck et al. 2010) has been used to test 
our results and calculate corrections to our abundances 
line by line. As shown in Fig. [51 the scattering corrections 
are negligible for red lines, but can approach 0.4 dex in 
the blue. Lines with the largest corrections are the res- 
onance lines of Mn I. Notice also that the corrections 
are sensitive to the atmospheric parameters, in particu- 
lar metallicity; the three stars in Carina that are shown 
were choosen because they have the largest corrections 
and they are the most metal-poor. 

We compared the elemental abundances between the 
programs MOOG and CALRAI (Spite 1967, Cayrel et 
al. 1991, Hill et al. 2012), and also the scattering correc- 
tions calculated with MOOG-SCAT to those determined 
with the program Turbospectrum (Alvarez & Plez 1998). 
The standard abundance results and scattering correc- 
tions were nearly identical for all lines between these 
codes, in particular, the scattering corrections were in 
agreement to within ~0.01 dex. Only the resonance lines 
of Mn I showed slightly larger differences in the scattering 
corrections (~0.04 dex); since resonance lines form over 
more atmospheric layers, we expect these lines are more 
sensitive to the details in the models and the line forma- 
tion calculations between the codes. Thus, we consider 
the consistency in the abundance results from the line 
formation codes MOOG and CALRAI, and the scatter- 
ing corrections between MOOG-SCAT and Turbospec- 
trum to be in excellent agreement. 

As a final note, the scattering corrections had no effect 
on our spectroscopic temperature or microturbulence de- 
terminations, nor the metallicities for [Fe I/H], since only 
Fe I lines with A > 5000 A were used in those steps. 

3.8.3. Gravity is not from Fe I/Fe II 

The analysis often stars in Carina by Koch et al. (2008) 
and five more stars by Shetrone et al. (2003) used the 
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Fig. 5. — A comparson of the line abundance corrections 

from MOOG-SCAT for one of the M68 standard star and the 

three most metal-poor stars in Carina. The y-label "Corr" = 
MOOG-SCAT - MOOG. 



Fe II lines to determine spectroscopic gravities from ion- 
ization equilibrium (where [Fe I/H] ~ [Fe II/H]). For 
four stars in common with Koch et al. , our gravities 
are smaller than theirs by log 0.5 dex. Since turbu- 
lent velocity and ionization equilibrium are correlated in 
RGBs, then the higher microturbulent values found by 
Koch et al. (see Section 13.5.21) are correlated with their 
higher spectroscopic gravities. 

There are sufficient Fe II lines in many of the stars 
in our analysis to examine the ionization equilibrium 
of iron and determine spectroscopic gravities, however 
we do not use this method in this analysis. For consis- 
tency with the FLAMES/GIRAFFE analysis, where a 
shorter wavelength interval meant fewer Fe II lines, then 
we adopt photometric gravities using the same methods 
as in Lemasle et al. (2012). On the other hand, we notice 
that the Fe I and Fe II abundances are in good agree- 
ment for our Carina targets, with a mean [Fe II/Fe I] = 
+0.02, and range of -0.24 to -^0.20, which is -lcr(Fe I). 
The stars in M68 have higher offsets, with a mean [Fe 
II/Fe I] = -1-0.35, which may reflect uncertainties in its 
distance or reddening. When the abundance of Fe II is 
found to be larger than Fe I in red giants, it can be due 
to the overionization of Fe I by the radiation field which 
is neglected under the assumption of LTE. Corrections 
for this effect can be as large as 0.3 dex (Mashonkina et 
al. 2010), which is similar to the offset in the M68 stars. 
Thus, we consider our Fe I and Fe II results to be in good 
agreement for all of our targets when physical gravities 
are adopted. 

Note, Fe II lines at all wavelengths were examined in 
this analysis, i.e., they were not restricted to A < 5000 
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TABLE 6 

Abundance Uncertainties for M68-6023 





Elem 




Alogg 




A[Fe/H] 


Added in 






+55 


-0.04 


+0.12 


-0.02 


Quadrature 


— 
- 


Fe I 


+0.09 


0.00 


—0.05 


0.01 


0.10 


- 


TT 

re 11 


—0.02 


—0.02 


—0.05 


—0.01 


0.06 


- 


OM 


+0.05 


—0.04 


0.00 


0.00 


0.06 


- 


Na I 


+0.08 


+0.01 


—0.03 


0.01 


0.09 




Mg I 


+0.08 


+0.02 


—0.03 


0.01 


0.09 


- 


Oa 1 


+0.08 


+0.01 


—0.04 


0.01 


0.09 




TT 

DC 11 


+0.02 


—0.01 


—0.05 


—0.01 


0.06 




Ti I 


+0.09 


+0.01 


—0.04 


0.01 


0.10 




Ti II 


+0.02 


—0.01 


—0.05 


—0.01 


0.06 




V I 


+0.06 


0.00 


-0.03 


0.01 


0.07 




Cr I 


+0.10 


+0.01 


-0.05 


0.01 


0.11 




Mn I 


+0.10 


0.00 


-0.07 


0.01 


0.12 




Co I 


+0.10 


+0.01 


-0.08 


0.01 


0.13 




Ni I 


+0.09 


+0.01 


-0.05 


0.01 


0.11 




Zn I 


+0.01 


0.00 


-0.02 


0.00 


0.02 




Sr II 


+0.04 


-0.01 


-0.08 


0.01 


0.09 




Y II 


+0.02 


-0.01 


-0.03 


-0.01 


0.04 




Zr II 


+0.02 


-0.02 


-0.02 


-0.01 


0.03 




Ba II 


+0.04 


-0.02 


-0.09 


-0.01 


0.10 




Eu II 


+0.03 


-0.02 


-0.09 


-0.01 


0.10 



-3 



-2.5 -2 
[Fel/H] 



-1.5 



-1 



Fig. 6. — Iron abundances ([Fel/H]) from high resolution spectral 
analyses versus those derived from the CaT calibration by Starken- 
burg et al. (2010). A line of equal motallicities is shown (dashed), 
±0.2 dex (dotted). 

A like Fe I. 

3.9. Evaluating the CaT Metallicity 

The initial inetallicities in our analysis are taken from 
a calibration of the near-IR Ca II triplet near 8500 A(= 
CaT). A direct comparison between the low resolution 
CaT metallicity index and high resolution iron abun- 
dances for large samples of RGB stars in the dwarf galax- 
ies Fornax and Sculptor was shown by iBattaglia et al.l 
((2OO81) to provide consistent abundances in the range 
—2.5 < [Fe/H] < —0.5. More recently, Starkenburg et 
al. (2010) found that the metallicities deviate strongly 
from the linear relationship for RGB stars with [Fe/H] 
< —2.0 and have developed a new calibration for metal- 
poor stars that also considers the [Ca/Fe] ratio. 

The CaT metallicities for our targets are listed in Ta- 
ble [5l and plotted against the [Fe I/H] values in Fig. [6l 
The high resolution [Fe I/H] abundances tend to be 
slightly larger than the metallicities from the CaT index, 
however, as Starkenburg et al. (2010) show, the ratio of 
[Ca/Fe] becomes more important in the calibration of low 
metallicity stars. The two most metal-poor stars in this 
analysis are found to have lower [Ca/Fe] ratios than the 
Galactic trend which was assumed in the Starkenburg et 
al. calibration. 

4. ABUNDANCE RESULTS 

4.1. Solar Abundances and Galactic Comparisons 

For compari son purposes, solar abundances from 
lAsplund et al.l (2009) are adopted. Galactic comparisons 
are from the compilation by Venn et al. (2004) , supple- 
mented with Si, V, Cr, and Nd from Fulbright (2000). 
The compilation was updated with thick disk stars from 
Reddy et al. (2006), La from the compilation by Roed- 
erer et al. (2010), and Cu and Zn from Mishenina et 



TABLE 7 

Abundance Uncertainties for CAR.-1087 



Elem 


AT.ff 


Alogg 




A [Fe/H] 


Added in 




+98 


-0.10 


-0.25 


-0.06 


Quadrature 


Fe I 


+0.10 


0.00 


+0.08 


+0.01 


0.13 


Fe II 


-0.03 


-0.02 


+0.08 


-0.01 


0.09 


Na I 


+0.13 


0.00 


+0.06 


+0.01 


0.15 


Mgl 


+0.05 


+0.02 


+0.04 


+0.01 


0.07 


Ca I 


+0.08 


+0.01 


+0.02 


+0.01 


0.09 


Sc II 


+0.03 


-0.02 


+0.06 


-0.01 


0.07 


Ti I 


+0.17 


+0.01 


+0.08 


+0.01 


0.19 


Ti II 


0.00 


-0.02 


+0.08 


-0.01 


0.08 


V I 


+0.13 


+0.02 


+0.04 


+0.01 


0.14 


Cr I 


+0.18 


+0.01 


+0.09 


+0.02 


0.20 


Mn I 


+0.18 


+0.02 


+0.10 


+0.02 


0.21 


Ni I 


+0.07 


+0.01 


+0.03 


+0.01 


0.08 


Ba II 


+0.08 


-0.04 


+0.02 


-0.01 


0.09 


Eu II 


+0.04 


-0.03 


+0.02 


-0.01 


0.06 



al. (20 02). All of these were rescaled to lAsplund et al.l 
(|2009() solar abundances from the Grevesse & Sauval 
(1998) values. The compilation of abundances of metal- 
poor stars in the Galactic halo and Local Group dwarf 
galaxies by Frebel (2010b, after eliminating the upper 
limit values) was also used here for comparison purposes. 
Previously published abundances fo r Car ina and M68 
were re-scaled to the lAsplund et al.l (|200i^ ) solar abun- 
dances, i.e., Shetrone et al. (2003) and Lee et al. (2005) 
from Grevesse fc Sauval ( 1989), and Koch et al. (2008) 
from lAsplund et aD ()2005D values. 

4.2. Abundance Error Estimates 

Errors in the abundances have been determined in two 
ways, with the maximum of these errors adopted for the 
analysis (Tables [6] to [8]) . Firstly, the error in the equiva- 
lent width measurements are conservatively estimated as 
AEW = EWrms + 10% X EW, where EW^ms is from the 
Cayrel formula (Section l3.2l) . This AEW was propagated 
through the abundance analysis to provide a a{EW) in 
each line abundance. This error is not necessarily sym- 
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TABLE 8 

Abundance Uncertainties for Car-484 



Elem 




Alog a 


A? 


A[Fe/H] 


Added in 




+66 


-0.1 


+0.07 


—0.02 


Quadrature 


Fe I 


+0 06 


—0.01 


—0 04 


0.00 


0.07 


Fe II 


—0.07 


—0.05 


—0.03 


0.00 


0.09 


O I 


0.00 


—0.05 


—0.01 


—0.01 


0.05 


Na I 


+0.06 


+0.01 


—0.01 


0.00 


0.06 


Mo- T 

IVig i 


+0.05 


+0.01 


—0.03 


0.00 


0.06 


Q: T* 
Ol i 


0.00 


—0.01 


0.00 


—0.01 


0.01 


Ca I 


+0.09 


+0.01 


—0.03 


0.00 


0.10 


Sc II 


—0.02 


—0.03 


—0.02 


0.00 


0.04 


T^: T 
ill 


+0.13 


—0.01 


—0.04 


0.00 


0. 14 


T^: TT 
i 1 ii 


—0.03 


—0.05 


—0.04 


0.00 


0.07 


V 1 


+0.14 


—0.01 


—0.01 


0.00 


0.14 


Cr I 


+0.14 


-0.01 


-0.05 


0.00 


0.15 


Mn I 


+0.10 


0.00 


-0.03 


0.00 


0.11 


Co I 


+0.08 


-0.02 


-0.03 


0.00 


0.08 


Ni I 


+0.04 


-0.01 


-0.03 


0.00 


0.05 


Cu I 


+0.08 


-0.02 


-0.04 


0.00 


0.08 


Zn I 


-0.04 


-0.02 


-0.03 


0.00 


0.05 


Y II 


0.00 


-0.04 


-0.03 


0.00 


0.05 


Ba II 


+0.01 


-0.06 


-0.06 


0.00 


0.08 


La II 


+0.01 


-0.06 


-0.01 


-0.01 


0.06 


Nd II 


0.00 


-0.05 


-0.03 


0.00 


0.06 


Eu II 


-0.01 


-0.05 


-0.01 


-0.01 


0.05 



Si I errors determined from the lines and EWs in Car-524. 



metric, and we adopt the average. For synthetic spectra, 
(j{EW) was estimated from the range of abundances for 
which a good fit of the observed hne profile could be 
achieved. The second method for estimating the error 
in the abundances is simply to calculate the standard 
deviation when more than five lines of an element are 
available, cf{X); when fewer lines are available, we set 
a{X) = cr(Fe I). The final error in [X/H] is adopted as 
the maximum of a{X) or a{EW). The corresponding 
error on the mean, S{X), is calculated as 

6iX) = a{X)/^ 

where Njc is the number of lines available of element 
X. To get the final error in [X/Fe], i.e., the ratio with 
iron, then the error on the mean ((5(FeI)) was added in 
quadrature. 

These line measurement errors do not take into account 
the errors due to the uncertainties in the stellar parame- 
ters. These errors are difficult to ascertain and combine 
properly since they are not independent and their rela- 
tionship(s) to one another are not well defined. We re- 
port these errors for three representative stars, Car-6023, 
Car-1087, and Car-484, in TablesEHH andEl These three 
stars represent the range in stellar parameters, resolu- 
tion, and SNR of all of our spectra. It can be seen that 
the stellar parameter errors are similar to or smaller than 
the line measurement errors (exceptions include Ca I, 
Ti I, Cr I, and Co I where the stellar parameter errors 
can be slightly larger, but still <0.15 dex). 

The stellar parameter errors reported for Car-6023 
can be applied to Car-6022 and Car-6024, while those 
for Car-1087 can be applied to Car-5070 and Car-7002. 
Those for Car-484 can be applied to the rest of the 
FLAMES/UVES sample, ahhough the lower SNR for 
Car-1013 means the errors for that target are larger by 
~2x (as seen in Table [5]). 



TABLE 9 

M68 Chemical Abundances 



Elem 


6022 




6023 


6024 




[X/Fe] ±cr ( 


;#) 


[X/Fe] ±<T (#) 


[X/Fe] ±<T (#) 


Fe I 


4.99 ±0.02 


(92) 


5.09 ±0.03 (92) 


4.86 ±0.03 (81) 


Fe II 


0.40 ±0.05 


(19) 


0.26 ±0.05 (19) 


0.39 ±0.06 (17) 


CH 


-0.82 ±0.20 


(S) 


-0.62 ±0.20 (S) 


-0.59 ±0.20 (S) 


O I 


< 1.36 


< 1.33 


< 1.54 


Na P 


0.07 ±0.10 


(S) 


0.27 ±0.10 (S) 


-0.03 ±0.10 (S) 


Mg I 


0.01 ±0.20 


(S) 


-0.09 ±0.20 (S) 


0.04 ±0.20 (S) 


Si I 


< 1.22 


< 1.15 


< 1.36 


Ca I 


0.48 ±0.04 


(17) 


0.39 ± 0.04 (16) 


0.48 ±0.05 (15) 


Sc II 


0.57 ±0.05 


(10) 


0.28 ± 0.10 (10) 


0.44 ±0.06 (9) 


Ti I 


0.11 ±0.05 


(23) 


0.12 ± 0.04 (21) 


0.30 ±0.07 (20) 


Ti II 


0.77 ±0.05 


(52) 


0.54 ± 0.05 (51) 


0.65 ±0.05 (46) 


V I 


-0.06 ±0.16 


(2) 


0.02 ± 0.17 (2) 


0.35 ±0.17 (2) 


Cr I 


-0.17 ±0.06 


(12) 


-0.25 ± 0.09 (12) 


-0.17 ±0.08 (9) 


1\/Tn 


n AS -1-0 07 


yi) 


n P.^ -1-0 ^'y (7\ 
— u.Di diu.iz yt J 


n 1 -1-0 1 {A\ 
— v.oi mu.iu 


Co I 


-0.27 ±0.09 


(6) 


-0.26 ±0.06 (7) 


0.10 ±0.10 (6) 


Ni I 


0.06 ±0.03 


(8) 


0.20 ±0.10 (8) 


0.04 ±0.14 (7) 


Zn I 


0.47 ±0.16 


(2) 


0.33 ±0.17 (2) 


0.43 ±0.24 (1) 


Sr II 


-0.08 ±0.16 


(2) 


-0.26 ±0.17 (2) 


-0.23 ±0.17 (2) 


Y II 


-0.07 ±0.23 


(1) 


-0.43 ±0.17 (2) 


-0.21 ±0.17 (2) 


Zr II 


0.20 ±0.23 


(1) 


-0.05 ±0.24 (1) 


0.00 ±0.24 (1) 


Ba II 


0.33 ±0.12 


(4) 


0.14 ±0.12 (4) 


0.36 ±0.12 (4) 


La II 


< 1.37 


< 1.36 


< 1.56 


Nd II 








1.27 ±0.24 (1) 


Eu II 






0.24 ±0.17 (2) 


0.47 ±0.17 (2) 


Note. 


— We calculate [X/Fe] = [X/H] - [Fe 


I/H]. For Fe I, the 



abundance listed is 12+log(Fc/H). 

^ A correction of A[Na/Fe] = —0.5 has been applied to account for 
NLTE effects on the Na D lines (Andrievsky et al. 2007). Similarly, 
we apply a correction of Alog (Mn/H) = +0.5 to the Mn I reso- 
nance lines to account for NLTE effects (Bergemann & Gehren 2008). 



4.3. M68 Comparison Stars 

The stellar parameters and chemical abundances for 
our three stars in 1VI68 are compared to the one star anal- 
ysed by Shetrone et al. (2003) and eight stars by Lee et 
al. (2005). The targets in this analysis are fainter by '-^ 1.4 
magnitudes, thus they are located further down the RGB 
and have slightly higher temperatures and gravities. Our 
microturbulence values are comparable or lower. 

The metallicities from Fe I for these three stars are in 
excellent agreement with one another, as expected for 
a single stellar population. The mean Fe I abundance 
is log(Fe I/H) = 5.04 ±0.10, where the error is the line 
to line abundance scatter in each star, added in quadra- 
ture, and divided by v^- This is similar to the mean 
[Fe I/H] abundance determined by Lee et al. (2005), 
where [Fe I/H] = 4.96 ±0.06 when they adopt the photo- 
metric gravities. Higher iron abundances are found with 
lower spectroscopic gravities; this is also seen by Shetrone 
et al. (2003) for one star. 

Most of the elemental abundances determined here are 
in excellent agreement with those from Lee et al. (2005) 
and Shetrone et al. (2003), and our Na I values are within 
the ranges of the others. To compare the [X/Fe] ratios 
for these elements with those from Lee et al. requires 
knowing that their [X/Fe] results use Fe I for ratios with 
neutral species and Fe II for ratios with ionized species 
and oxygen In our analysis, we use only the Fe I lines 
to indicate the metallicity, thus we adjust their [X/Fe II] 
abundances by adding their [Fe II/Fe I] ratio, = +0.40. 

The most significant differences are that our Mg and 



12 



Venn et al. 



TABLE 10 

Carina Chemical Abundances 



TABLE 11 

Carina Chemical Abundances continued 



Eiem 




484 


















705 






769 






1013 






[X/Fe] ±a 


(#) 


[X/Fe] ±<j 


(#) 


[X/Fe] ±a 


(#) 




[X/Fe] ±a 


(#) 


[X/Fe] ±o- 


(#) 


[X/Fe] ±a 


(#) 


Fe I 


5.97 


±0.02 


(^^A') 


5 75 _(_n no 


(121) 




("1 261 


Fe I 


6 


15 


±0.02 


(127) 


5.82 


±0.02 


(125) 


6.20 


±0.03 


(125) 


Fe II 


0.07 


±0.04 




n 09 +n 03 


(15) 


n 1 2 +n 05 


(18) 


Fc II 





02 


±0.06 


(1&) 


0.06 


±0.05 


(15) 


—0.24 


±0.06 


Ol) 


O I 


0.40 


±0.21 


(W 

y^i 


n 41 -|_n oi 


(1) 


< 0.09 




I 


<0 


12 












<0.08 






Na I 


—0.25 


±0.21 




< —0.26 


< —0.72 




Na I 


— 


22 


±0.25 


(1) 


0.04 


±0.25 


(1) 


—0.25 


±0.37 


(1) 


Mp- 1 


0.19 


±0.15 




n 27 +n 1 5 


(2) 


_n Kn _l_n ic 
yj.ijyj 1 w ■ 1 1 1 


(2) 


Met T 





13 


±0.18 


(2) 


0.42 


±0.18 


(2) 


0.09 


±0.26 


(2) 


Si I 






«n _|_n 1 r: 


C?) 


n 1 Q _|_n OQ 


(^) 


Si I 





23 


±0.25 


^^ ) 


0.45 


±0.18 


(2) 


0.06 


±0.26 


(9) 


Ca I 


0.16 


±0.04 


(20) 


n 1 n _|_n 04 


yin) 


_n 1 7 _|_n 04 


("191 
yii)) 


Ca I 





11 


±0.04 


f19~) 
yLVj 


0.24 


±0.06 


(2^) 


—0.06 


±0.05 


(2n) 


Sc II 


0.06 


±0.04 


C12'1 
yLZ,) 




yiv) 




(8) 


Sc II 


-0 


03 


±0.05 


(11) 


0.15 


±0.06 


(11) 


-0.24 


±0.07 




Ti I 


0.23 


±0.04 


(25) 


_n no -Un 0"^ 


C?"?) 
yzz) 


_n 49 -1-0 04 


yLO) 


Ti I 





02 


±0.05 


(2K\ 


0.27 


±0.05 


(2X\ 


0.27 


±0.06 


(2A) 


Ti II 


0.18 


±0.05 




n HQ _|_n nfi 


(13) 


_n on _|_n ne; 


(12) 


Ti II 


-0 


01 


±0.06 


(11) 


0.13 


±0.06 


(11) 


0.11 


±0.09 


(^ X) 


V I 


0.01 


±0.03 


(10) 


0.21 ±0.15 


(2) 


-0.67 ±0.23 


(1) 


V I 





02 


±0.06 


(10) 


0.41 


±0.11 


(5) 


0.09 


±0.08 


(11) 


CrI 


0.00 


±0.12 


(4) 


-0.09 ±0.10 


(4) 


-0.20 ±0.12 


(4) 


Cr I 


-0 


23 


±0.12 


(4) 


0.08 


±0.12 


(4) 


0.25 


±0.21 


(3) 


Mn I 


-0.27 


±0.05 


(7) 


-0.34 ±0.10 


(4) 


-0.51 ±0.10 


(6) 


Mn I 


-0 


30 


±0.07 


(7) 


0.02 


±0.11 


(5) 


-0.51 


±0.17 


(5) 


Co I 


-0.13 


±0.21 


(1) 


-0.27 ±0.21 


(1) 


< -0.86 


Co I 


-0 


26 


±0.25 


(1) 


0.02 


±0.25 


(1) 


-0.01 


±0.26 


(2) 


Ni I 


-0.07 


±0.04 


(15) 


-0.08 ±0.04 


(14) 


-0.46 ±0.04 


(14) 


Ni I 


-0 


04 


±0.04 


(17) 


0.03 


±0.04 


(13) 


-0.10 


±0.06 


(17) 


Cu I 


-0.47 


±0.21 


(1) 


< -0.85 




< -1.60 




Cu I 


-0 


80 


±0.25 


(1) 


-0.55 


±0.25 


(1) 








Zn I 


0.16 


±0.21 


(1) 


-0.17 ±0.21 


(1) 


-0.83 ±0.23 


(1) 


Zn I 


-0 


05 


±0.25 


(1) 






-0.46 


±0.37 


(1) 


Y II 


-0.32 


±0.12 


(3) 


-0.46 ±0.12 


(3) 


-1.37 ±0.16 


(2) 


Y II 


-0 


63 


±0.18 


(2) 


-0.27 


±0.14 


(3) 


-0.32 


±0.21 


(3) 


Ba II 


0.16 


±0.26 


(1) 


0.09 ±0.21 


(1) 


-0.64 ±0.16 


(S) 


Ba II 


-0 


61 


±0.14 


(3) 


-0.17 


±0.14 


(3) 


-0.15 


±0.26 


(S) 


La II 


0.48 


±0.12 


(3) 


0.55 ±0.15 


(2) 


< 0.11 


La II 


<0 


15 












<0.14 






Nd II 


0.63 


±0.15 


(2) 


0.52 ±0.15 


(2) 


-0.56 ±0.23 


(1) 


Nd II 





26 


±0.18 


(2) 


0.27 


±0.18 


(2) 


0.11 


±0.26 


(2) 


Eu II 


0.47 


±0.21 


(1) 


< 0.63 


< 0.26 


Eu II 


<0 


31 




<0.63 




0.33 


±0.37 


(1) 



Note. — We calculate [X/Fe] 
abundance listed is 12+log(Fe/H). 



[X/H] - [Fe I/H]. For Fe L 



the Note. — We calculate [X/Fc] 
abundance listed is 12+log(Fe/H). 



[X/H] - [Fe I/H]. For Fe I, the 



IVIn abundances are lower than Lee et al. 's (2005), and 
the Shetrone et al. 's (2003) values for Zn, Y, Ba, and 
Eu are lower (see Figures in the next Sections). For the 
heavy elements examined by Shetrone et al. , we track 
the differences to the model atmosphere parameters, e.g., 
ionized species that are more sensitive to gravity, and 
Zn abundances that are ~2x more sensitive to temper- 
ature than the other elements. For IVIn, we note that 
Lee et al. examined only a single, very weak line (IVIn I 
6021.8), whereas ours are from several lines that are in 
good agreement with one another, thus we consider our 
result more reliable. For IVIg, our abundances have been 
determined from four spectral lines — three lines from 
spectrum synthesis due to severe blending and difficulties 
in the continuum placements (two near 3830 A, and one 
at A5172), and a fourth line at A5528 which is sufficiently 
unblended to be analysed from EW measurements (note 
that a fifth line is available in the FLAMES/UVES spec- 
tra at A5711 but is too weak to be measured in any of the 
JVIagellan/lVIIKE spectra). The IVIg abundances are con- 
sistent between these four lines and from star to star in 
1VI68, however we assign a larger uncertainty to account 
for the difficulties in setting the continuum in the three 
line syntheses. The average [IVIg/Fe] is slightly lower than 
found in the other analyses (Lee et al. 2005, Shetrone et 
al. 2003), than the Galactic trends at that metallicity, 
and than found in the other a-elements. While the offset 
is within the estimated uncertainties, we note that our 
[Na/Fe] are slightly higher than the Galactic trends — 
it is possible these stars exhibit some deep mixing (e.g., 
Gratton et al. 2010). 

New elemental abundances not previously determined 
in JVI68 (Sr II and Zr II) are in good agreement with the 
Galactic distributions. Finally, our upper limits for Si I, 
O I, and La II are in agreement with the abundances de- 
termined by Lee et al. (2005) and Shetrone et al. (2003). 



Elem 



TABLE 12 
Carina Chemical Abundances 



1087 
[X/Fe] ±<7 (#) 



5070 
[X/Fe] ±<7 (#) 



7002 
[X/Fe] ±<T (#) 



Fe I 


4.69 ±0.06 (38) 


5.35 ±0.05 (55) 


4.64 ±0.05 


(38) 


Fe II 


0.01 ±0.07 (7) 


-0.14 ±0.10 (11) 


0.20 ±0.14 


(6) 


CH 


< -1.02 (S) 


< -0.98 (S) 


< -1.07 (S) 


O I 


< 1.55 


< 0.90 


< 1.61 






Na P 


-0.82 ±0.24 (S) 


-1.18 ±0.38 (S) 


-0.45 ±0 


23 


(S) 


Mg I 


0.48 ±0.34 (1) 


-0.36 ±0.38 (1) 


0.22 ±0 


23 


(2) 


Si I 


< 1.70 


< 1.01 


< 1.76 




Ca I 


-0.03 ±0.20 (3) 


-0.01 ±0.10 (9) 


0.18 ±0 


10 


(6) 


Sc II 


-0.01 ±0.15 (5) 


-0.20 ±0.24 (8) 


0.10 ±0 


15 


(4) 


Ti I 


0.65 ±0.17 (4) 


-0.08 ±0.15 (8) 


0.11 ±0 


18 


(3) 


Ti II 


0.29 ±0.07 (24) 


0.04 ±0.07 (30) 


0.39 ±0 


08 


(20) 


V I 


0.72 ±0.34 (1) 


-0.26 ±0.38 (1) 








Cr I 


-0.82 ±0.17 (4) 


-0.47 ±0.08 (7) 


-0.34 ±0 


16 


(4) 


Mn I 


-0.72 ±0.24 (2) 


-0.81 ±0.18 (3) 


-0.33 ±0 


23 


(2) 


Co I 


< 1.11 


< 0.31 


-0.26 ±0 


23 


(2) 


Ni I 


0.36 ±0.34 (1) 


0.18 ±0.13 (7) 


0.35 ±0 


18 


(3) 


Cu I 


< 0.37 


< -0.43 


< 0.39 




Zn I 


< 0.66 


< 0.01 


< 0.77 






Sr II 


-1.56 ±0.24 (S) 




-1.21 ±0 


32 


(S) 


Y II 


< -0.05 


< -0.69 


< -0.12 




Zr II 


< 0.24 


< -0.41 


< 0.36 






Ba II 


-1.07 ±0.34 (S) 


-1.12 ±0.27 (S) 


-0.98 ±0 


18 


(S) 


La II 


< 1.77 


< 1.15 


< 1.75 




Nd II 


< 0.99 


< 0.34 


< 1.09 






Eu II 


< 1.85 


0.13 0.38 (1) 


< 1.90 







Note. — We calculate [X/Fe] = [X/H] - [Fc I/H]. For Fe I, the 
abundance listed is 12+log(Fe/H). 

^ A correction of A[Na/Fe] = —0.6 has been applied to account for 
NLTE effects on the Na D lines (Andricvsky ct al. 2007). Similarly, 
we apply a correction of Alog (Mn/H) = +0.5 to the Mn I reso- 
nance lines to account for NLTE effects (Bergemann & Gehren 2008) . 
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Wavelength 

Fig. 7. — Spectrum synthesis of the CH 4320 A feature in 
M68-6023. The chemical composition determined in this paper 
was used for the syntheses, along with three carbon abundances, 
[C/Fe] = —0.4, —0.6, —0.8. Broadening parameters were initially 
taken from the resolution of the Magellan spectra (Gaussian 
FWHM ~ 0.15) and checked against the best fit well defined 
iron-features. 

Thus, our M68 resuKs appear to be in good agree- 
ment with other abundance determinations for this chis- 
ter, and Galactic cluster and field stars in general. We 
conclude that our stellar parameter determinations and 
model atmospheres analyses can therefore be used reli- 
ably for elemental abundances in the Carina targets. 

4.4. Carina Stars 

The chemical abundances of the stars in Carina are 
presented in Tables [10] to [121 and discussed element by 
element in the following Sections. 



^ 
u 



CEMP 



-4 



-3 



-2 
[Fc/Hl 



-1 



CEMP 




12 3 4 

Log (L/Lsun) 

Fig. 8. — Carbon abundances in Galactic stars (small grey dots; 
Frebel 2010b, Gratton et al. 2000), the three stars in M68 (empty 
black squares), and upper limits on the three metal poor Carina 
stars {upper panel). CEMP stars are usually defined as having 
[C/Fe] > 1.0, although extra mixing at the RGB tip can lead to 
a luminosity dependence as noted in the lower panel. The low 
carbon abundances in the Carina stars are consistent with those 
of luminous RGB stars in globular clusters (data from Gratton et 
al. 2000), which are not similar to CEMP stars. 



4.4.1. Carbon 

Carbon forms during the helium burning phases, 
whether the hydrostatic helium core burning phases in 
massive stars, or helium shell burning phases in AGB 
stars (Woosley & Weaver 1995). The chemical evolution 
of carbon is further complicated by CNO cycling, where 
carbon is reduced when exposed to hot protons and the 
CN-cycle runs to equilibrium values. This evolution in 
carbon can be seen in RGB stars in globular clusters, 
where the [C/H] abundance and ^^C/ C ratio are both 
reduced as stars ascend the RGB due to internal (self) 
mixing (Gratton et al. 2004) . 

Carbon abundances in this analysis were determined 
from spectrum syntheses of portions of the A^H - X^A 
CH G-band near A4320 in the Magellan spectra. The 
best fit was taken as the carbon abundance from the CH 
line list by Brown et al. (1987) and Carbon et al. (1982), 
with acceptable fits as the abundance uncertainty. A 
(^i2^Qi3 j-atio of six was adopted (the typical value for 
bright RGB stars that have undergone the first dredge 
up) ; this ratio made only a small difference such that in- 
creasing the value to 50 caused Alog(C/H) — 0.06. Sim- 
ilarly, reducing the oxygen abundance had only a small 
effect (since C can be locked into the CO molecule); we 
adopted [0/Fe] = 4-0.4, but reducing this by 3x lowered 
the carbon abundance by only < 0.05 dex. The spectrum 
of M68-6023 around the CH 4320 A feature is shown in 
Fig. [71 with spectrum syntheses for three carbon abun- 
dances. 

The carbon abundances for the M68 stars are in good 
agreement with typical red giants in globular clusters 



above the RGB bump (e.g.. Smith et al. 2005, Gratton 

et al. 2000), where [C/Fe] 0.4 at the RGB bump and 

slowly decreases to ~ —1.0 dex at the RGB tip due to 
mixing of CNO-cycled gas. This trend can be seen in 
Fig. [8] (lower panel, data from Gratton et al. 2000). The 
carbon abundances compiled by Frebel (2010b) suggests 
a much larger scatter in carbon in metal-poor stars than 
seen in clusters (Fig. [8l upper panel), however Frebel's 
compilation includes field stars of all evolutionary stages, 
including carbon enhanced metal-poor (CEMP) stars. 
CEMPs have [C/Fe] > 1.0 and are thought to be both 
stars that have undergone mass transfer in a binary sys- 
tem with an AGB companion (CEMP-s, e.g.. Beers & 
Christlieb 2005) and stars with high primordial carbon 
(e.g., from h igh m ass rotating stars, Meynet et al. 2006). 
lAoki et all (20071) suggested that mixing on the RGB 
and extra mixing at the tip of the RGB will lower the 
surface carbon abundance, such that the definition of 
a carbon enhanced metal-poor (CEMP) star should de- 
pend on luminosity. This luminosity dependent range 
is shown in Fig. [51 (lower panel), and clearly shows that 
stars in M68 and the Carina dSph are not CEMP stars 
by any definition. 

The mean carbon abundance in M68 is <[C/Fe]> — 
— 0.7± 0.1, in agreement with typical upper RGB stars. 
Upper limits on the carbon abundances in the metal- 
poor Carina stars are even lower ([C/Fe] < —1.0); this 
value is typical of upper RGB stars in globular clusters, 
but on the low side of the values found for the bright, 
metal-poor field stars in the Galaxy, as seen in Fig. [HI 



14 



Venn et al. 



4.4.2. Alpha Elements 



The alpha-elements (O, Mg, Si, Ca, Ti) are built from 
multiple captures of He nuclei (a particle) during various 
stages in the evolution of massive stars (> 8 Mq, e.g., 
carbon burning, neon burning, complete and incomplete 
Si burning), and dispersed during SN II events. Thus, the 
[a/Fe] ratio in a star is a way to trace the relative contri- 
butions from SN II to SN la products in the ISM when 
it formed (the nucleosynthesis of iron is discussed in the 
next Section). Although Ti is not a true alpha-element, 
the dominant isotope '^^Ti forms through explosive Si- 
burning and the a-rich freeze oufl thus it behaves like 
an alpha-element (Woosley & Weaver 1995). 

Our [a/Fe] results are plotted in Figs. l9l and fTOl Koch 
et al. (2008) determined a element abundances for 4 stars 
in common with ours, and a solid line connects those 
results in the Figures. In general, the a-element abun- 
dances in Carina tend to be lower than the Galactic dis- 
tributions at all intermediate metallicities, i.e., —2.2 < 
[Fe/H] < —1.3. In the most metal-poor stars at [Fe/H] = 
—2.9, most of the a elements are in good agreement with 
the Galactic distributions. Two stars stand out in their 
[a/Fe] distributions, Car-612 and Car-5070, which have 
very low ratios of [Mg,Ca,Ti/Fe] for their metallicities. 
These two stars will be discussed further in the following 
Sections. 

Oxygen: Oxygen abundances in this analysis are solely 
from the [O I] 6300 line. Oxygen could only be measured 
in the higher metallicity stars in the sample where the 
line is stronger. Excellent agreement was found for one 
star in common with Koch et al. (2008; Car-484, where 
they found [0/Fe] =4-0.39 and we find =-f0.48). 

Silicon: Silicon abundances are only available for the 
higher metallicity stars. Line measurements of the Si I 
lines at 5684 A and 6155 A are similar to those in Koch et 
al. (2008), however our oscillator strength for one line is 
lower by 0.4 dex. Only upper limits were available from 
the Magellan spectra taken of the metal-poor stars. 

Magnesium: Mg I abundances were calculated from 
the A5528 line in all of the Carina stars, combined with 
A5172 and A5183 when <200 mA, and A5711 in the 
FLAMES /UVES spectra (it was too weak to measure in 
the Magellan spectra). The abundances between these 
lines are in good agreement. We note that the Carina 
Mg abundances are in good agreement with those from 
Koch et al. (2008) and Shetrone et al. (2003), see Fig.\M 

Calcium: Many lines of Ca I were available in all of the 
program stars, over a wide range of wavelengths, and of 
appropriate line strengths. Results are in good agree- 
ment with Koch et al. (2008), see FigfTOl 

Titanium: Titanium abundances were determined 
from a range of Ti I and Ti II lines across the entire spec- 
tral region. The abundances between the two species are 
in good agreement, typically within 1 cr of each other, 
despite the possibility that titanium can be overionized 
by the radiation field resulting in lower Ti I abundances 
when LTE is assumed. Letarte et al. (2010) suggested 
that better agreement is found when [Ti I/Fe I] is com- 

® The a-rich freeze out occurs when a strong shock travels 
through infaUing matter at the inner most shell of a core-collapse 
supernova. After being shocked, material is heated to > 5 billion 
degrees, so that nuclei decompose into neutrons and tightly bound 
a particles. 
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Fig. 9. — The LTE oxygen and silicon abundances of stars in 
Carina and M68, compared to the Galactic distribution. Results 
for Carina are from this paper (filled circles), Koch et al. (2008, 
empty circles), and Shetrone et al. (2003, filled triangles); for 
four stars in common between this analysis and Koch et al. , the 
abundance results are connected by a grey line. M68 results are 
shown for one star from Shetrone et al. (2003, grey triangle), the 
mean of 7 RGBs from Lee et al. (2005, grey square). Usually 
our M68 results are shown by empty black squares, but for O 
and Si only upper limits are available, shown as grey arrows. A 
representative error bar is shown (A [Fe/H] = ±0.1, A[X/Fe] = 
±0.15), although the actual errors per star from this analysis 
are also plotted per point. Small grey points show the Galactic 
distributions summarized by Venn et al. (2004) and Frebel (2010b). 

pared to [Ti II/Fe II] in their sample of Fornax RGB 
stars, due to similarities in the NLTE effects and pos- 
sibly temperature scale variations. Examination of our 
M68 comparison stars are in agreement with their find- 
ing, such that the agreement improves when Ti II is com- 
pared with Fe II, however this had no effect on our Carina 
abundances. Thus, in this paper, the abundances from 
both species are averaged together. We note that our 
Carina Ti abundances are in good agreement with those 
from Koch et al. (2008) and Shetrone et al. (2003); see 

Fig.iini 

4.4.3. Sodium 

Sodium is produced primarily during the carbon burn- 
ing stages, but as metallicity increases then a sufficient 
amount of neutron rich material also allows sodium to 
be produced through the Ne-Na cycle during H-burning 
(Woosley & Weaver 1995). Thus, the chemical evolution 
of Na initially follows the a elements, but then deviates 
from this (rising) when AGB stars contribute. Further- 
more, Herwig et al. (2004) have shown that metal-poor 
AGB stars from 2 to 5 Mq will overproduce sodium by 
factors of 10 to 100. This Na overproduction has multiple 
sites dependent on the AGB mass. 

In the FLAMES/UVES spectra, sodium abundances 
were determined from a single Na I line at 5688 A; this 
was because the resonance Na D lines at 5895 and 5889 
A were too strong (~ 300 mA), and the line at 5682 
A (used by Koch et al. 2008) was too weak. Only the 
most metal-poor stars in Carina have Na D lines that 
are weak enough to be used in this analysis, but since 
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Fig. 10. — The LTE magnesium, calcium and titanium abun- 
dances of stars in Carina and M68, compared to the Galactic 
distribution. Symbols are the same as in Fig. [9l 



those stars were observed with Magellan/MIKE spectra, 
their lower SNR spectra required spectrum syntheses for 
the line abundance determinations. The Fe I 5615 and 
5634 lines were examined to check the broadening param- 
eters (Gaussian broadening with FWHM = 0.25 A was 
adopted, which is the effective resolution of the Magel- 
lan MIKE spectra). When the Na D lines could not be 
fit simultaneously, then an average of the best fit abun- 
dances for each line was adopted, and the range used to 
es timate the uncerta i nty. 

(jAndrievskv et all I2007D have examined the NLTE 
corrections for the Na D resonance lines in metal-poor 
RGB stars. Using their Table 2, combined with their 
Fig. 6 showing the effect of the correction on the line 
equivalent widths, then we estimated a NLTE correction 
for our abundances. For the M68 stars, the higher gravi- 
ties and strong line strengths imply corrections of ~ —0.5 
dex. For the three metal-poor Carina stars (Car-1087, 
Car-5070, and Car-7002) then the atmospheric parame- 
ters and line strengths suggest corrections of ^ —0.6 dex. 
These corrections are included in Fig. [TT] Of course, no 
corrections are applied to the other Carina objects since 
we did not use the Na D lines for those. 

For Car-612, we only determined an upper limit, EW 
< 30 mA for Na I A5688, which corresponds to a very 
low upper limit of [Na/Fe] = —0.8. This result is compa- 
rable to that from Koch et al. (2008) who measured an 
EW (= 24 mA) for the same weak line and determined 
[Na/Fe] — —0.5. Our lower upper limit is due to small 
differences in the atmospheric parameters and adopted 
solar Na abundance. 

4.4.4. Iron peak Elements 



^ 

S5 



CAR »Oi 

M68 □ □ A " 



7002 
1087 



'■I 



-3 2 -1 

[Fel/II] 

Fig. 11. — Sodium abundance ratios of stars in Carina and M68, 
compared to the Galactic distribution. For the metal-poor stars in 
Carina and M68, NLTE corrections have been applied to the Na 
D line abundances; the symbols are the same as in Fig. (9] For the 
Galactic data, o nly a limited amount of available data are shown; 
the results from JAndrievskv et al.ll2007l ) and Gehren et al. (2006) 
arc included since these two studies considered NLTE effects on the 
NaD lines, and results from Reddy et al. (2003, 2006) are shown 
since they only used the subordinate Na I lines. 

In the early Universe, the iron-peak elements (Sc to 
Zn: 23 < Z < 30) are exclusively synthesized during 
Type II supernovae explosions by explosive oxygen and 
neon burning, and complete and incomplete explosive Si 
burning. Abundances of these elements show a strong 
odd-even effect (odd nuclei have lower abundances than 
the even nuclei) . Their yields depend on the mass of the 
progenitors (e.g., Woosley & Weaver 1995), the mass cut 
adopted (mass expelled relative to the mass that falls 
back onto the remnant, e.g., Nakamura et al. 1999) and 
SNe II explosive energies (e.g., Umeda & Nomoto 2005). 
Only at later times (~ 1 Gyr, e.g., Maoz et al. 2010), 
when lower mass stars reach the end of their life time, 
do SNe la become a significant, possibly the dominant, 
contributor to the total iron-group inventory. The on- 
set of SNe la in the chemical evolution of the Galaxy is 
observed as a knee in the [a/Fe] vs [Fe/H] near [Fe/H] 
-1.0 (e.g., McWilliam et al. 1995). At lower metal- 
licities, SNe la [X/Fe] yields could vary (e.g., Kobayashi 
& Nomoto 2009). We have measured several of the Fe- 
peak elements in the Carina stars to compare with the 
Galactic trends. 

Scandium: Several (4-13) lines of Sc II were measured 
and HFS corrections applied. The distribution in the Sc 
abundances at intermediate metallicities is larger than 
most of the other iron peak elements (see Fig. [12] and 
[T5|) . although this includes Car-612 which has several 
chemical peculiarities and stars analysed by Shetrone et 
al. (2003). Other than those objects, the Sc abundances 
are similar to the Galactic distribution. 

Vanadium: For most of the intermediate metallicity 
stars, the V abundance is well determined from ~10 V I 
lines and HFS corrections were applied. However, in the 
metal-poor stars and Car-612, it is determined from only 
1-2 lines. 

Manganese: The Mn I abundances were determined 
from several lines (4-7) in the intermediate metallicity 
stars, and these line abundances showed good agreement 
with one another. Generally, [Mn/Fe] in Carina follows 
the Galactic distribution above [Fe/H] = — 2. In the 
metal-poor Carina stars and M68 stars, a different set of 
spectral lines was used (because of the wavelength cover- 
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Fig. 12. — The odd-Z LTE abundance ratios for scandium, 
vanadium, and manganese in Carina. HFS corrections have been 
applied, as well as a NLTE correction to the resonance lines of 
Mn I. Symbols are the same as in Fig.|9]virith the exception of the 
Galactic standards for [Mn/Fe] which are taken from Sobeck et 
al. (2006) and Cayrel et al. (2004). 

age and resolution of the Magellan/MIKE spectra). The 
blue resonance lines at A4030, 4033, and 4034, as well as 
several additional subordinate lines were measured, with 
only A4823 in common between the MIKE/Magellan and 
FLAMES/UVES datasets. The Mn results from the res- 
onance lines were lower than from A4823 and the other 
subordinate lines in the M68 standard stars, most likely 
due to non-LTE effects. Bergemann & Gehren (2008) 
have calculated non-LTE corrections for Mn I lines and 
do find large corrections for the blue resonance lines that 
are metallicity dependent. They do not give corrections 
for red giants, but for red main sequence stars the cor- 
rections can be as large as -1-0.5 dex at [Fe/H] = —3. 
These corrections are similar to the offsets that we find 
between resonance and subordinate lines, therefore we 
apply a correction of Alog(Mn) — +0.5 dex to the Mn I 
resonance line abundances — this correction should be 
checked with proper modelling of these metal-poor red 
giant atmospheres. Correcting the resonance line abun- 
dances improves the mean [Mn/Fe] results, but the mean 
values in Car-5070 and Car-1087 remain lower than for 
similar stars in the Galactic halo. 

Chromium: Cr abundances were determined from 4-7 
lines of Cr I, which showed good agreement from line to 
line. A recent calculation of the NLTE effects on the Cr I 
lines by Bergemann & Cescutti (2010) showed that the 
corrections are small (< 0.1 dex), but that Cr I/Cr II 
ionization equilibrium and the solar [Cr/Fe] ratio is re- 
gained for metal poor stars, rather than the downward 
trend seen in Fig. [T31 Our Carina data generally follow 
the Galactic trend, with the exception of very low abun- 
dances in Car-5070, Car-1087, Car- 705, and Car-612. 



Fig. 13. — LTE abundance ratios for chromium, cobolt, and 
nickel in Carina. The low abundances for Car-612, and low Cr in 
Car-5070, are exceptional, like in Fig. 1121 Symbols are the same 
as in Fig. |9] 

Cobalt: Co was determined from only 1-2 lines in 
this analysis. The [Co/Fe] ratios lie slightly below the 
Galactic abundances, although they are not well con- 
strained from so few lines. The only exception is the 
very low upper lim it to [Co/Fe] for Car-612 (see Fig. lTS)) . 
iBergemann et alT] ([2010) calculated NLTE corrections 
for Co I and Co II; they suggest the NLTE corrections 
depend on metallicity and can become as large as -1-0.6 

to +0.8 dex at [Fe/H] 3.0. This would not affect the 

very low Co upper limit found for Car-612. 

Nickel: Several Ni I lines (13-17) were measured in 
the intermediate metallicity stars over a range of wave- 
lengths, showing that [Ni/Fe] is in agreement with the 
Galactic distribution (Fig. [T3| . Only 1-3 Ni I lines were 
available in the most metal-poor stars though, thus the 
uncertainties are large because a(Fe I) was adopted for 
the error estimates. Only the Ni result for Car-612 stands 
out, but this result is robust, determined from 14 Ni I 
lines with a small error in the mean. Combined with 
the low Na upper limit, then Car-612 fits the Na-Ni re- 
lationship observed for some a-poor stars in the Galaxy 
by Nissen & Schuster (1997, 2010). 

4.4.5. Cu & Zn 

While contiguous on the Periodic Table, the main nu- 
cleosynthetic sites of copper and zinc are difficult to as- 
certain. In massive stars, Cu and Zn form during com- 
plete Si-burning, the a-rich freeze out, and even the weak 
s-process (e.g., Timmes et al. 1995). In SN la models, the 
yields of Cu and Zn are sensitive to the neutron excess 
and thus metallicity (Matteucci et al. 1993, Travaglio et 
al. 2005, Kobayashi & Nomoto 2009). The contributions 
to the production of Cu and Zn in AGB stars are uncer- 
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Fig. 14. — The LTE copper and zinc abundances in the Carina 
stars. The Cu and Zn abundances in Car-612 are remarkably low. 
Symbols are the same as in Fig. (9] 

tain, though recent calculations suggest small amounts 
of Cu and even smaller amounts of Zn can be produced 
in more massive (5 Mq) AGB stars (Karakas et al. 2008, 
Karakas 2010). Precise estimates of the AGB yields can 
also depend on uncertain parameters such as the mass 
loss law and number of dredge up episodes (Travaglio et 
al 2004). 

Copper: The Cu I line at 5105 A was observed in some 
of our target stars. Upper limits were calculated adopt- 
ing an EW of 30 mA for the FLAMES /UVES spectra and 
50 mA for the most metal-poor Carina stars. Hyperfine 
structure corrections were applied. Generally, copper fol- 
lows the Galactic trend, and though we have very few 
actual measurements, the upper limit for Car-612 is very 
low and provides a very strong constraint. 

Zinc: The Zn I line at A4810 was observed in most 
of our sample, or used to determine the upper limits 
to the Zn abundance (adopting <40 mA); see Fig [T4l 
The zinc abundances at intermediate metallicities were 
slightly lower than the Galactic distribution, especially 
Car-612. 

4.4.6. Neutron Capture Elements 

Neutron-capture elements (Sr to U, 38 < Z < 92) orig- 
inate in the rapid neutron capture process (r-process) 
that occurs during explosive nucleosynthesis in SNe II. 
The r-process has a main source in 8-10 Mq stars that 
form elements with Z>50, and a weak source in more 
massive progenitors (>20 Mq) that contribute primar- 
ily to the lighter neutron-capture elements (e.g., Sr, Y, 
Zr, Travaglio et al. 2004). A more detailed examina- 
tion of the r-process shows that a more accurate rep- 
resentation is in terms of a high entropy wind model 
(Farouqi et al. 2009, Roederer et al. 2010); nucleosynthe- 
sis in low entropy winds proceeds primarily through a 
chared-particle (a-) process, whereas a neutron-capture 
component (the classical weak and main r-processes) oc- 
cur in the high entropy winds. Thermal pulsing in AGB 
stars also contributes to these elements through the slow 



neutron capture process (s-process). With intermediate 
masses (2-4 Mq), AGB stars have longer lifetimes than 
the sites of the r-process, and therefore in simple chem- 
ical evolution models AGB stars do not produce any of 
the heavy elements in the most metal-poor stars, but 
dominate the formation of some of these elements at later 
times and higher metallicities. In the Galaxy, this is seen 
as the rise in the s-process w hich begins ne ar [Fe/H] = 
-2.5 (McWilliam et al. 1998, iBurris I [20Q0l Francois et 
al. 2007); this is in contrast to the knee in [a/Fe] due to 
contributions from SNe la near [Fe/H] = —1.0. In the 
Sun, the r-process contributes (11%, 15%, 25%, 28%, 
53%, and 97%o) of ( Sr, Ba, La, Y, Nd, and Eu) respec- 
tively (|Burris II2000D . Thus Eu is critically important as 
a nearly pure r-process indicator. 

Yttrium: [Y/Fe] in the intermediate metallicity stars 
was determined from Y II 4883, 5087, and 5200, which 
gave very consistent abundances from line to line and 
star to star. The Y II 4900 line could not be used be- 
cause it is severely blended in all of the spectra. The 
[Y/Fe] ratios are similar to those found by Shetrone et 
al. (2003), i.e., slightly below the Galactic values at in- 
termediate metallicities, see Fig. [151 One exception is 
Car-612 which has an extremely low [Y/Fe] result, and 
another is Car-5070 which has a low upper limit value. 
The low [Y/Fe] in Car-612 can be seen directly in com- 
parison with Car-705, where Fig. ll6l which shows that the 
Fe I line strengths are similar but the Y II 5087 line is 
much weaker in Car-612. Upper limits were determined 
in the most metal-poor stars using 40 mA as the EW 
upper limit for the Y II A4883 line. 

Barium: Five lines of Ba II were analysed (at 4554, 
4934, 5853, 6141, and 6496 A), which yielded consistent 
abundances from line to line and star to star, although 
they were not all observed in one star. Spectrum syn- 
theses were used to confirm the abundances when the 
SNR was low. Hyperfine structure corrections for three 
lines (5853, 6141, and 6496 A) are negligible (< 0.02). 
NLTE corrections are also negligible (< 0.03, Short & 
Hauschlidt 2006), with the exception of the A4554 reso- 
nance line, however we have chosen to not correct that 
line since the estimated correc tion is not large (< 0.15 
dex. Short & Hauschhdt 2006, lAndrievskv et al. 1120091) 
and the LTE abundance from that line is consistent with 
the other Ba II line results. Most stars in Carina have 
the same [Ba/Fe] distribution as stars in the Galaxy, with 
the exception of three stars, Car-612, Car-5070, and Car- 
705; see Fig. [HI 

Lanthanum: La was determined from three lines of 
La II at A6320, 6390, 6774 (additional lines at A4333, 
5301, 5303 could not be used due to the SNR of the spec- 
tra). Negligible HFS (< 0.02) were calculated for two La 
fines (A6320 and 6390). [La/Fe] follows the Galactic dis- 
tribution for the few stars where we could measure it (see 
Fig. [15]). 

Strontium: [Sr/Fe] was determined from two Sr II lines 
at very blue wavelengths (4077 and 4215 A) reached only 
by our Magellan/MIKE spectra. The SNR is poor in 
that region, thus we adopted the results from spectrum 
syntheses for those lines in the three Carina stellar spec- 
tra, see Fig. [TH] Nevertheless, the interpretation of the 
synthetic results remains difficult; the A4077 line is the 
stronger of the Sr II resonance lines, and yet the detec- 
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Fig. 15. — The LTE abundances for the heavy elements Y, Ba, 
and La in Carina and M68, compared to the Galactic distribution. 
The Y abundance in Car-612 is remarkably low, but verified by a 
comparison of the spectra of Car-612 and Car- 705, e.g., in Fig. 1161 
The Galactic star abundances for La are from the critically exam- 
ined compilation by Roederer et al. (2010). These ratios are again 
quite low for Car-612, and Ba in Car-5070. Symbols are the same 
as in Fig. |9] 
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Fig. 16. — Spectra around the Y II 5087 A lines in Car-612 and 
Car-705. Since the atmospheric parameters and metallicities of 
these two stars are very similar, then this figure shows yttrium 
is truly weaker in Car-612. Spectrum synthesis in this region 
confirms the different Y II abundances derived from the EW 
analyses. The lower [Ni/Fe] abundances in Car-612 can also be 
seen directly by comparing the 5084 A lines. 

tion of the A4215 line is more clear in Car-7002, and 
possibly Car-1087. The best estimates for the Sr abun- 
dances in these two stars are listed in Table [T^ A Sr 
upper limit for Car-5070 could not be determined due to 
noise near both the A4215 and 4077 hncs. NLTE cor- 
rections were considered, but not applied to our results. 
Short & Hauschildt (2006) estimate abundance cor rec- 
tions of —0.07 dex, whereas (|Andrievskv et aLll2011[ ) es- 
timate corrections of ~ 0.0 and +0.1 dex for A4077 and 
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Fig. 17. — LTE abundances for the heavy elements Sr, Nd, and 
Eu in Carina and M68, compared to the Galactic distribution. 
Symbols the same as in Fig. O 

4215, respectively. We chose to not apply either cor- 
rection since these estimates are small and in opposite 
directions. 
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Fig. 18.— Spectrum syntheses of the Sr II 4077 and 4215 
lines in the three most metal-poor Carina stars, for [Sr/Fe] = 
—0.5, —1.0, —1.5, and — 2.0. The low SNR makes accurate Sr abun- 
dances very difficult. 

Zirconium: Upper limits for Zr were determined from 
the Zr II A4208 line in the Magellan spectra. These upper 
limits are quite high, and do not add the discussion on 
the heavy element abundances. 

Neodymium: Nd was determined from two Nd II lines 
(A5319 and 5249) for all stars with FLAMES/UVES 
spectra, and these lines were used to calculate the up- 
per limits from the Magellan spectra. Nd is consistent 
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Fig. 19. — Neutron capture ratios with Eu of stars in Carina 
and M68, compared to the Galactic distribution. The ratio with 
Eu permits an assessment of the s-process contributions relative 
to the r-process contents in the Sun (Burris et al. 2000, dashed 
line). The rise in the s-process can be seen in the Carina stars and 
the Galactic distributions. Upper and lower limits are shown for 
some stars. Symbols the same as in Fig.fO] 



with the Galactic distribution, though it is very low in 
Car-612 (Fig.HT]). 

Europium: The europium abundance in the Sun is 
nearly entirely due to the r-process, therefore Eu is iden- 
tified as an important indicator in any stellar analysis to 
establish the ratio of r-process to s-process contributions 
of the neutron capture elements. In this analysis, three 
Eu II lines were examined (A4129, 4205, & 6645). Un- 
fortunately, all three lines were not observed in any one 
star; the Magellan spectra have too low resolution to de- 
tect the weak 6645 A fine, and the FLAMES/UVES data 
do not cover the bluer wavelengths. Nevertheless, the Eu 
results are consistent with the Galactic distribution (see 
Fig I17p . Hyperfine structure and isotopic splitting cor- 
rections are negligible (<0.05 dex). 

4.4.7. s-process to r-process Ratios: [X/Eu] 

To assess the relative amounts of s-process to r- 
process contributions amongst the neutron capture el- 
ements, then abundance ratios relative to Eu are exam- 
ined. In Fig. [m [Y/Eu] and [Ba/Eu] of stars in Ca- 
rina and the Galaxy are shown and compared to the so- 
lar r-process fractions from Burris et al. (2000; also see 
lArlandini et~al] 119991 In both [Y/Eu] and [Ba/Eu] (as 
well as [La/Eu] and [Nd/Eu], not shown), the gradual 
rise in these ratios can be seen with increasing metallic- 
ity. McWilliam et al. (1998) showed that the rise in the 
s-process starts at [Fe/H] — —2.5 in the Galaxy. In Ca- 
rina, this rise appears to begin at a higher metallicity, 
[Fe/H] > —2.0. This is most likely due to the metallicity 
dependence of the AGB yields, as modelled by Travaglio 
et al. (2004, also Pignatari et al. 2010), where fewer iron 



seed nuclei in a high neutron density wind can collect 
more neutrons, thus underproducing the first s-process 
peak elemental, and overproducing the second and/or 
third s-process peak elements. The slightly lower [Y/Eu] 
and slightly higher [Ba/Eu] ratios in Carina are consis- 
tent with metal-poor AGB stars contributions. 

As a final note, the high Ba abundance in M68 is some- 
what surprising. This was also found by Lee et al. (2005). 
Since [Ba/Eu] is high, but [Eu/Fe] is normal, this sug- 
gests M68 has had an unusual chemical evolution, possi- 
bly having been enriched (or self-enriched) by metal-poor 
AGB stars. 

5. DISCUSSION 

Detailed chemical abundances for up to 23 elements 
in nine stars in the Carina dSph have been presented in 
this paper. Previously, only five stars in this dSph had 
detailed abundance determinations for so many elements 
(Shetrone et al. 2003), while another ten had [a/Fe] de- 
terminations from high resolution spectroscopy (Koch et 
al. 2008). In a companion paper by Lemasle et al. (2012), 
the VLT FLAMES/GIRAFFE spectra for 36 additional 
RGB stars are presented, but due to the limited wave- 
length coverage of those spectra, only four elements are 
analysed in detail (Mg, Ca, Ba, and Fe). For the first 
time, Lemasle et al. (2012) calculated the ages of all 
stars with detailed chemical abundances in the Carina 
dSph, including those in this paper, from isochrone fit- 
ting. While the uncertainties in the ages can be quite 
large for any one star, the differential ages should be bet- 
ter, and any discussion regarding ages is therefore limited 
to two age bins, an old population and an intermediate- 
aged population. These ages (two age bins) are adopted 
in this discussion. 

In the previous Section, our results were compared to 
stars in the Galaxy, where some chemical peculiarities 
were noted. In this Section, we will compare our abun- 
dance results to those of stars in other dwarf galaxies. 
Three dSphs have been chosen for comparison: Sculptor, 
Fornax, and Sextans. The datasets for the three dSphs 
are taken fro m Shetr one et al. (2001, 2003), Geisler et 
al. f2005). iAoki et al.l ([200a) . Letarte et al. (2006, 2010), 
Tafelmeyer et al. (2010), and Frebe l et al. (2010c), all 
scaled to the lAsplund et al.l (|2OO90 solar abundances. 
These three dSphs were chosen since large datasets of de- 
tailed abundances from high resolution spectroscopy are 
available, and the analyses of those spectra used similar 
methods and model atmospheres to our analysis. We do 
not include resuhs from Kirby et al. (2009, 2010, 2011) 
in these comparisons to avoid systematic differences due 
to the lower resolution of their spectra and differences 
in the analysis methods. Also, we do not include higher 
mass dwarfs (e.g., the LMC, SMC, or Sgr) because the 
metallicity range of the stars in those galaxies do not 
overlap well with the stars in the Carina dSph. 

We also compare several elemental abundances of the 
most metal-poor stars in Carina to stars in the ultra 
faint dwarf galaxies (UFDs) and the metal-poor stars 
in the Draco dSph. Datasets for the UFDs are from 

The s-process peaks are defined by the neutron magic numbers 
for filling nuclear shells (e.g., N=50, 82, 126 are full nuclear shells). 
A full nuclear shell lowers the cross section for further neutron 
captures, thus elements collect at these neutron numbers defining 
the first, second, and third s-process peaks. 
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the compilation by Frebel (2010b), and for Draco are 
from Cohen et al. (2009), Fulbright (2004), and Shetrone 
et al. (2001). These comparisons are valuable because 
Carina has a tiny dynamical mass within its half light 
radius (M1/2), e.g., M1/2 (Carina) = 6.1 ±2.3 x 10^ M© 
(Walker et al. 2009), which is a factor of 2, 4, and 9 times 
smaller than Sculptor, Sextans, and Fornax, respectively, 
but is comparable to those of the more massive UFDs. As 
an example, Bootes I has M1/2 = 5.9 ±3.7 x 10^ M©, 
which is essentially the same as Carina, although the 
masses of the UFDs continue to be revised downwards, 
e.g., the mass of Bootes I has been revised downwards 
by a factor of 14 due to the lower velocity dispersions 
found by Koposov et al. (2011). Of course. Carina is not 
considered an UFD because of its luminosity - Carina is 
8x brighter than Bootes I (Walker et al. 2009). 

5.1. The Metallicity Distribution of the High Resolution 
Spectroscopic Sample 

As shown in Fig. 4 of Lemasle et al. (2012), the 
stars that have been analysed with high resolution spec- 
troscopy are biased towards metallicities near [Fe/H] 
= —1.5 ± 0.3. This is the same as the mean metal- 
licity and metallicity range predic ted for the domi- 
nant intermediate-aged population bv lBono et al.l ()2010f ) 
from their CMD analysis; however this is not sci- 
entifically significant for two reasons. Firstly, more 
than half of the high resolution sample comes from 
the FLAMES/GIRAFFE spectral analysis presented by 
Lemasle et al. (2012), and as discussed in that paper, the 
low SNR of those spectra meant that the weaker lines of 
the more metal-poor stars were harder to detect. This 
effectively removed those stars from the analysis. Sec- 
ondly, although the V magnitude at the tip of the RGB 
may be slightly brighter with higher metallicities, this is 
only a small effect and maximizing the fiber placement 
was a more significant concern. Also, while it is true that 
the three stars targeted for Magellan/MIKE observations 
are located in the outer fields of Carina and are the most 
metal-poor stars sampled, this does not imply a popu- 
lation gradient in Carina as those stars were purposely 
selected for their low metallicities. Koch et al. (2006) 
found no significant difference in the mean metallicities 
of RGB stars in the inner and outer fields of Carina. 
Thus, the results presented in this paper are not used to 
constrain the distributions in location or metallicity of 
the stars in Carina. 

5.2. Dispersions in [a/Fe] 

In Fig. [201 we show [Mg/Fe] and [Ca/Fe] for Carina, 
compared to the three other dSphs and the Galaxy, over 
the full metallicity range examined. Looking at [Ca/Fe] 
alone suggests that the chemical evolution of Carina has 
been similar to Sculptor. In Sculptor, a noticeable down- 
ward trend in these abundances begins at [Fe/H] ~ —1.8, 
whereas in Galactic stars this occurs at a higher metallic- 
ity. This knee in the [a/Fe] ratios is usually interpreted 
as the onset of contributions to [Fe/H] from SN la, with 
the shift in the knee to lower metallicities in dwarf galax- 
ies attributed to their slower chemical evolution (e.g., 
Lanfranchi et al. 2008, Kirby et al. 2011). In Carina, it is 
not clear where or if there is a knee. 

The more remarkable result seen in Fig. [20] is 
the dispersion in [Mg/Fe] (A[Mg/Fe]) in Carina, 



which is observed in both the FLAMES/UVES and 
FLAMES/GIRAFFE data, and is much larger than the 
dispersion in [Ca/Fe] (A[Ca/Fe]). A calculation of the 
intrinsic spreadf[3 suggests that this difference is real 
and signficant when considering the whole dataset, the 
data from Lemasle et al. (2012) only, our nine UVES 
targets alone, or the data at a specific intermediate 
metallicity (such as [Fe/H]= -1.2 ± 0.1), i.e., 0.4 < 
[N, ( [Mg/Fe] )-Ni( [Ca/Fe])] < 0.7. Only for [Fe/H] < 
—2.0 are the intrinsic spreads between [Mg/Fe] and 
[Ca/Fe] similar. 

Differences in the dispersions of [Mg/Fe] and [Ca/Fe] 
may be partially due to differences in their nucleosyn- 
thetic sites, but also may be the result of inhomogeneous 
mixing of the interstellar gas and therefore poor statisti- 
cal sampling of the SN contributions when forming stars. 
In terms of nucleosynthesis, SN yields of Ca and espe- 
cially Mg depend on the progenitor mass (e.g., Woosley 
& Weaver 1995, Iwamoto et al. 1999); for example, Mg 
forms in hydrostatic core C and O burning, whereas Ca 
has contributions from the a-rich freeze out and explo- 
sive Si burning during the SN II explosion Differences 
in the SN la models can also lead to differences in the 
[Si-Ca/Fe] ratios, e.g., the central density, metalhcity, ig- 
nition source, fiame speed, and even type of SN la model 
can play a role (e.g., Maeda et al. 2010, Ropke et al. 2006, 
Iwamoto et al. 1999). In terms of inhomogeneous mix- 
ing, models by Revaz & Jablonka (2011, also Revaz et 
al. 2009) predict a large spread in the [Mg/Fe] ratios in 
low mass dwarf galaxies due to longer gas cooling times 
and subsequently longer mixing timescales for the inter- 
stellar medium. As an example, the hot gas from SNe II 
can be subject to buoyant forces requiring up to 2 Gyr 
to cool and mix through a low mass galaxy. This hot gas 
can also quench star formation. Therefore, it is possible 
to find models of low mass dwarf galaxies that predict 
both a high dispersion in element ratios and an episodic 
star formation history, like Carina. 

Surprisingly, the [Mg/Ca] ratios in Carina do not show 
a larger dispersion then the other dSphs in Fig. [501 This 
may imply that inhomogeneous mixing plays the domi- 
nant role (over differences in nucleosynthetic sites). The 
one star with the extremely low [Mg/Ca] ratio is Car- 
743, analysed by Lemasle et al. (2012) from their lowest 
SNR spectrum; the formal uncertainty in that one result 
is ~ 2x larger than the representative errorbar shown, 
and should be considered with caution. The one star 
with the highest [Mg/Ca] ratio is Car-1087, analysed in 
this paper; with a metallicity of [Fe/H] = —2.9 and high 
[Mg/Ca] ratio then this star is very similar to the unusual 
star Draco-119 (discussed below). 

5.2.1. Metal-poor Stars and the UFDs 

In Fig.[2Tl we show only the metal-poor tail in the Ca- 
rina abundances and compare the [Mg/Fe] and [Ca/Fe] 
ratios to those of metal-poor stars in five UFD galaxies 
and the Draco dSph. Most of the UFD stars have high 
values of [Mg/Fe] and [Ca/Fe], like the Galactic halo, 
whereas Carina and Draco show a range extending to 
very low values. 

We have estimated the intrinsic spread (N^) from the formula 
N? = - (cr)^, where N is the range in [X/Fe] and (ct) the average 
error in [X/Fe]. 
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Fig. 20. — Mg and Ca abundances of stars in Carina (red), 
Sculptor (green), Fornax (cyan), and Sextans (magenta). Symbols 
for Carina: solid circles are from this paper, squares are from 
Lemasle et al. (2012), triangles are from Shetrone et al. (2003), and 
asterisks are from Koch et al. (2008). Symbols for Fornax: circles 
are from Tafelmeyer et al. (2010), triangles are from Shetrone et 
al. (2003), squares are from Letarte et al. (2010). Symbols for 
Sculptor: large/small squares are UVES/GIRAFFE FLAMES 
data from Hill et al. (2012), solid circles are from Tafelmeyer 
et al. (2010), asterisk are from Giesler et al. (2005), triangles 
are from Shetrone et al. (2003), and the star is fro m Frebel et 
al. (20 10c). Symbols for Sextans: squares are from lAoki et al.l 
l|2009l ). triangles are from Shetrone et al. (2001), and circles 
are from Tafelmeyer et al. (2010). Tafelmeyer et al. (2010) note 
that the very low [Ca/Fe] ratio reported for the most metal-poor 
dSph star (at [Fe/H] ~ -4) is likely due to NLTE effects in the 
formation of the strong resonance Ca I 4227A line, a line that 
was not used throughout the rest of their analysis, and therefore 
should be regarded with caution. Representative error bars of 
A [Fe/H] = ±0.1 and A[X/Fe] = ±0.15 are shown. 

Some stars in all of these systems show very high 
[Mg/Ca] values (> 0.5) e.g., two stars in the Hercules 
UFD (Koch et al.2008), one in Bootes I (Feltzing et 
al. 2009), one in Draco (Fulbright 2004, Shetrone et 
al. 2001), and one in Carina (this paper). In the UFDs, 
it has been proposed that high [Mg/Ca] may be due to 
the chemical enrichment of its interstellar gas by as few 
as one SN II explosion from a massive progenitor (e.g., 
a >35 Mq star). In this scenario, a unique chemical 
signature can be imprinted onto the gas that is used to 
form stars, while the rest of the gas is expelled to quench 
further star formation (e.g., Koch et al 2008b, Frebel et 
al. 2010a, Simon et al. 2010). This scenario is aided by 
peculiar neutron capture ratios, e.g., in the Hercules stars 
the Ba II lines are not detected (Koch et al. 2008). Simi- 
larly, the Ba II (and Sr II) lines in the high [Mg/Ca] star 
in Draco (Draco-119) are also not detected (Fulbright et 
al. 2004 suggested that this star was enriched by a SN II 
in the mass range of 20 — 25 Mq , and is lacking in contri- 
butions from higher mass progenitors). However, Ba II 
(and Sr II) lines are detected in the high [Mg/Ca] star in 
Carina. Those lines are weak, leading to very low abun- 



FlG. 21. — [Mg/Fe] and [Ca/Fe] in the most metal-poor stars 
in the Galaxy, Carina, Draco, and UFD galaxies. The data 
for Carina are the same as in Fig. 1201 while for Draco, the 
small/large squares are from Shetrone et al. (2001)/Cohen et 
al. (2009). The large empty triangle is Draco-119, the metal-poor 
star found by Shetrone et al. (2001) and reanalysed with higher 
SNR spectra by Fulbright (2004). The UFDs are presented as 
follows: solid/empty circles are for Bootes I from Feltzing et 
al. (2009)/Norris et al. (2010), solid triangles are for Ursa Major 
II from Frebel et al. (2010a), solid diamonds are for Com Ber 
from Frebel et al. (2010a), plus signs are for Hercules from Koch 
et al. (2008b), and the asterisk for one star in Leo IV by Simon 
et al. (2010). The dashed line in the [Mg/Ca] panel represents 
stars with the more extreme enhancements (3 in UFDs, 2 in dSphs). 

dances of [Ba/Fe] and [Sr/Fe], but if the source of these 
elements is the same as in the other dSphs, then Carina 
appears to have been able to retain a bit more of those 
early (high mass progenitor) SN II enrichments. 

5.2.2. Measurement Errors? 

One question worth considering is whether the large 
dispersion in [Mg/Fe] is simply due to measurement er- 
rors. For example, there are only 1-2 Mg I lines avail- 
able in each star, whereas there are 3-21 lines of Ca I 
(depending on the metallicity) in this analysis. Sim- 
ilarly, due to the limited wavelength coverage of the 
FLAMES/GIRAFFE spectra, there are <13 Ca I hues 
measured by Lemasle et al. (2012), and only one Mg I 
line at A5528. A larger number of spectral lines reduces 
the random measurement errors in an elemental abun- 
dance (by \/N). This is reflected in the smaller average 
error in [Ca/Fe] than [Mg/Fe] in Table 10 of Lemasle et 
al. (2012). In this paper, even if the measurement error 
from fitting a line is small, as in a high SNR spectrum, 
the analysis of a single line means that the adopted er- 
ror is the standard deviation in the Fe I line abundance 
(cr(Fe); see Section H?^ . These (conservative) errors are 
about as large as the [Mg/Ca] dispersion, but smaller 
than the [Mg/Fe] dispersion, therefore the large disper- 
sion in [Mg/Fe] does not appear to be due to measure- 
ment errors alone. 
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Another consideration is that variations in the [Mg/Fe] 
and [Mg/Ca] ratios have been found in a number of 
metal-poor stars in the ultra faint dwarf galaxies, and 
attributed to incomplete mixing and poor sampling of 
the full mass function, as shown in Fig. [5TJ Only the 
largest outliers are examined in terms of pollution by a 
single SN II, and similarly large outliers are found in our 
[Mg/Ca] and [Mg/Fe] data as well. The consistency be- 
tween these analyses suggests that the signatures in the 
largest outliers are not due only to measurement errors. 

5.3. Age Considerations 

The Carina dSph is an interesting galaxy partially be- 
cause of its low mass and partially because of its episodic 
star formation history. Lemasle et al. (2012) were the 
first to attempt to interpret the chemical evolution of this 
dwarf galaxy after separating the stars into two popula- 
tions, old (> 10 Gyr) and intermediate-aged (< 6 Gyr); 
stars with ages between 6 and 10 Gyr were examined by 
Lemasle et al. 2012, but we do not include those in this 
discussion. The old population shows a large range in 
[Fe/H] and [Mg/Fe], whereas the intermediate-aged (lA) 
population appears to have a very small range in [Fe/H] 
and [Mg/Fe]. Lemasle et al. (2012) were also struck by 
the significant overlap in [Fe/H] between these popula- 
tions, which would imply that the lA stars formed from 
gas that was more metal-poor than the end-point metal- 
licity of the old population. Combined with the higher 
mean [Mg/Fe] abundance of the lA population, Lemasle 
et al. (2012) suggested that the second epoch of star for- 
mation in Carina may have occurred after the late accre- 
tion of metal-poor, a-element rich gas. 

The overlap in the metallicities of the old and lA pop- 
ulation can be seen in Figs. [H] to [Mj We note that 
the overlap is heavily weighted by the age assigned to 
Car-612, which is the chemically peculiar star discussed 
in Section 4. Differences in the specific chemistry of 
stars are known to affect their isochrone ages (Dotter et 
al. 2007), and thus the age assignment to this star should 
be considered with caution. The overlap is also heavily 
weighted by Car-524, which has a large uncertainty in its 
age assignment. If the ages for these two stars are ne- 
glected, then there is a sharp transition at —1.4 < [Fe/H] 
< —1.6 between the age groups, i.e., equivalent to the 
measurement errors in metallicity. Thus, the overlap in 
metallicity between the two age groups is not sufficiently 
clear to indicate an infall of metal-poor gas to form the 
second generation of stars. 

5.3.1. [a/Fe] Ratios Between the Age Groups 

An offset in the mean [Mg/Fe] abundance in the lA 
population relative to the old population can clearly be 
seen in Fig. [521 A Kolmogorov-Smirnov (K-S) test shows 
that the distribution of the [Mg/Fe] values for the old 
population is broader and with lower mean values, i.e., 
i [Mg/Fe] /^oid = 0.05 ±0.32 and i[Mg7Fe]iiA = 0.22 ±0.15, 
both with moderately high probabilities of having "nor- 
mal" distributions once the outlier in the lA dataset (at 
[Mg/Fe] = -0.95) is removed. This was also seen by 
Lemasle et al. (2012) for both [Mg/Fe] and [Ca/Fe], al- 
though the offset in [Ca/Fe] is not as large, which we 
confirm with a K-S test. We also note that these off- 
sets in the [Mg/Fe] and [Ca/Fe] ratios with age are still 
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Fig. 22. — [o/Fe] ratios of stars in the Carina dSph separated 
into age groups. Symbols the same as in Fig. 1201 though now 
red is used for the old population and blue for the lA population; 
ages are from Lemasle et al. (2012). Galactic comparison stars are 
shown as small grey dots (including data from Venn et al. 2004, 
Frebel et al. 2010b, and Reddy et al. 2003, 2006). Representative 
errorbars are shown based on the moan error in the Carina dataset. 

present when only the high resolution data are examined 
(this paper, Koch et al. 2008, and Shetrone et al. 2003), 
although they are not as clear. The same signature is also 
hinted at in the [Ti/Fe] ratios (from the high resolution 
data) , but there is insufficient data in the lA population 
for a meaningful K-S test. Since the large disperions 
in the [a/Fe] ratios discussed above (larger than seen in 
most dSph galaxies) are an indication of inhomogeneous 
mixing, the simplest explanation for the offset in [a/Fe] 
between the old and lA population is that the second 
epoch of star formation occurred in a-enriched gas. The 
small range in the [a/Fe] and [Fe/H] ratios suggests this 
gas was well mixed. 

5.3.2. The Rise m the s-Process 

Clearly the old and lA stellar populations in Carina 
show different chemical signatures, and it is interesting to 
investigate how the two populations might be connected. 
While we have discussed the a-elements, possibly the 
most valuable abundance ratios for studying the chemical 
evolution of a dwarf galaxy are of the neutron capture 
elements. 

In Fig. [231 it is possible to see that AGB stars have 
contributed to both age groups in Carina, but not un- 
tfi [Fe/H] > —2. This can be seen by the [Ba/Eu] ra- 
tios, where the lower dashed line represents the r-process 
contributions (in the Sun), such that stars with this ra- 
tio have been enriched by SNe II products only. Above 
[Fe/H] = —2, the [Ba/Eu] ratios slowly increase from the 
low r-process value through contributions to Ba from the 
s-process (the rise in the s-process). This is seen in the 
Galaxy, Fornax, and Sculptor as well, although the rise 
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Fig. 23. — The neutron capture abundance ratios of stars in 
the Carina dSph separated into age groups. Data references and 
symbols for Carina and the Galaxy are the same as in Figs. 1221 
Data for the other dwarf galaxies are noted only by colour 
(cyan for Fornax, gr een for Sculptor, magenta for Sextans) with 
references as in Figs. [20l 
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Fig. 24.— Abundance ratios of [Na/Fe], [Mn/Fe], and [Cr/Fe] 
for stars in the Carina dSph separated into age grou ps. E)ata 
references and symbols are the same as in Fig. 1221 w ith the 
excep tion of Na where Galactic data are from ifAndricvsk v et al.l 
I20071) . Gehren et al. (2006), Edvardsson et al. (1993), and Reddy 
et al. (2003, 2006), and also Mn where the Galactic data are from 
Sobeck et al. (2006) and Cayrel et al. (2004). 



begins at higher metalhcities in those systems. From the 
[Ba/Fe] ratios, it is also possible to see how the AGB 
contributions begin to dominate the nucleosynthesis of 
Ba, e.g., the dispersion in [Ba/Fe] at low metalhcities in 
the Galaxy is interpreted as the stochastic sampling of 
the (small amounts of) Ba from SNe II products early 
on, but this scatter lessens when the dominant contribu- 
tions from AGB stars come at later times. The flatness 
of the relationship in [Ba/Fe] at higher metalhcities in 
the Galaxy implies that the timescale and yields of Ba 
from the AGB stars and Fe from SNe la are similar, and 
the small scatter implies that the gas is well mixed. In 
Carina, [Ba/Fe] in the lA population may have a large 
scatter. If the large scatter is real, then this would sug- 
gest that the AGB contributions were not well mixed 
in the ISM at any age, which would be an interesting 
contrast to the uniformity in the [a/Fe] ratios in the I A 
population. However, several of these data points are 
from Lemasle etal (2012) where the measurements of a 
single Ba II line were very difficult due to the low SNR 
in that spectral region. The data from our high resolu- 
tion analysis is far more uniform above [Fe/H] — —2, as 
shown in Fig. [TSl 

Below [Fe/H] = —2 it is difficult to ascertain if there 
are any stars with AGB contributions. There is one 
star (Car-5070) with [Ba/Eu] below the r-process ratio, 
implying that this star shows no signs of AGB prod- 
ucts; however Carina was poorly mixed at early times, 
and metal-poor AGB stars could produce very little Ba. 
The high [Ba/Y] ratios in stars over [Fe/H] = — 2 show 
that most of the s-process elements came from metal- 
poor AGB stars, where the first s-process peak (Y) was 



bypassed in favour of the second s-process peak (Ba; 
as described further in Section I4.4.7P — at even lower 
metalhcities, both can be bypassed for the third s-process 
peak. Therefore we cannot ascertain the precise metal- 
licity when AGB stars began to contribute in Carina, but 
certainly there are contributions from metal-poor AGB 
stars in the stars with [Fe/H] > —2, and a hint of the 
rise in the s-process between -2.0 < [Fe/H] < —1.6. 

Further evidence for the rise in the s-process can be 
seen in the evoluti on in the [Na/Fe]; see F ig. [Ml Cayrel 
et al. (2004) and (jAndrievskv et al.l I2007D showed that 
the [Na/Mg] and [Na/Fe] ratios are flat and low in metal- 
poor stars in the Galaxy, suggesting that Na is produced 
with the a-elements in SNe II at a (low) fixed ratio. The 
Na abundance rises when AGB stars begin to contribute 
(see Section [4.4.3^ . We see a very similar trend in the 
[Na/Fe] ratios in Carina; in fact the rise in [Na/Fe] is 
very similar to the rise in the s-process elements. One 
significant difference though is that the initial [Na/Fe] 
value appears to be much lower in Carina than in the 
Galaxy. The final [Na/Fe] ratios may be similar in Carina 
as in the Galaxy (and Sculptor), and we see no offset in 
the Na abundances between the old and I A populations. 
A dispersion in Na is not clear in this small sample. 

It is interesting that the evolution of [Mn/Fe] and 
[Cr/Fe] are also similar to Na in Carina. The abundances 
of these elements are much lower than in the Galaxy at 
low metalhcities, but then rise above [Fe/H] = —2.0 to 
abundances that are similar to stars in the Galaxy with 
the same metallicity. There is also no difference in these 
elements between the old and lA populations. Like Na, 
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once AGB stars and SNe la begin to produce iron-group 
elements, those contributions will dominate over the ini- 
tial (low) SNe II yields. The timescale for AGB contri- 
butions is thought to be equal to or shorter than that of 
SNe la 1 Gyr, Maoz et al. 2010), thus it is likely that 
both AGB stars and SNe la begin to contribute to the 
ISM in Carina near [Fe/H] = -2.0. The SNe la contri- 
butions appear to be well mixed in the ISM, and there is 
no offset between the old and lA stars. 

5.3.3. Early Chemical Evolution by SNe II 

The earliest stages of chemical evolution in Carina can 
be examined from the element ratios of the stars below 
[Fe/H] = —2. As discussed above, these stars appear 
to have been enriched only in SNe II products. While 
there is a large range in the [a/Fe] and [Mg/Ca] ratios 
(when the whole dataset is examined, including previ- 
ously published results) , the ratios of Mn, Cr, and Na 
all start out very low but then rise to the Galactic values 
at intermediate metallicities. 

In addition, the [Sr/Fe] values are amongst the lowest 
of all stars analysed in the Galaxy, other dSphs, and 
the UFDs (see Fig. [25|. The [Sr/Fe] values for Car- 
1087 and Car-7002 are similar to the lowest values found 
for metal-poor stars in Draco, Bootes I, and Com Ber. 
Tafelmeyer et al. (2010) suggested that dwarf galaxies 
may have a lower floor in the [Sr/Fe] ratio (~ —1-2) 
when their mass is equal to or lower than Draco. The low 
[Sr/Fe] ratios in Carina lead to very low [Sr/Ba] ratios 
as well. [Sr/Ba] ratios are quite interesting because they 
tend to be high in metal-poor stars in the Galactic halo, 
which has been interpreted as evidence for an extra nu- 
cleosynthetic source in the formation of these elements in 
massive stars (Travaglio et al. 2004, Ishimaru et al. 2005, 
Qian & Wasserburg 2008, Farouqi et al. 2009, Pignatari 
et aL2010). 

The absence of a [Sr/Ba] enhancement at low metallic- 
ities in Carina suggests a lack of the excesses seen in the 
Galactic stars. It is not clear from which stellar mass 
range these excesses arise: Travaglio et al. (2004) sug- 
gest the main r-process occurs in 8-10 Mq stars, and 
that the excess must arise from the more massive SNe II, 
whereas Farouqi et al. (2009) suggest the excess occurs in 
SN II with lower entropy winds where a neutron-rich, a- 
freezeout can occur. The specific mass or energy range in 
Farouqi's models is not clear. If we follow Travaglio's sug- 
gestions, then the excess Sr may form in more massive, 
or higher energy SNe II (i.e., hypernovae), and there- 
fore those SNe II seem to be missing in Carina — either 
those stars did not form or their ejecta was driven out of 
CaringFl. 

Mn and Cr also form in hypernovae as the decay prod- 
ucts of complete and incomplete Si-burning (Umeda & 
Nomoto 2002, Nomoto et al. 2011), thus a lack of hy- 
pernovae might also explain the very low initial values 

Low upper limits on the [Sr/Fe] and [Ba/Fe] are also reported 
for Draco-119, however these cannot be used together to constrain 
[Sr/Ba], As an exercise, if the upper limits are taken as the actual 
values, then [Ba/Fe] = —2.6 and [Sr/Ba] = +0.1, which places this 
star squarely amongst the Galactic distribution. Thus, even the 
lack of a detection of the Ba II and Sr II spectra lines in Draco-119 
does not provide a strong constraint on the one shot hypernova 
model. The low values could be consistent with inhomogeneous 
mixing and/or gas driven out by supernova winds. 
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Fig. 25. — [Sr/Fe] and [Sr/Ba] in low metallicity stars in the 
Galaxy, classical dSphs, and UFDs. The [Sr/Ba] ratio is enhanced 
in metal-poor stars in the Galaxy, which is not always s een in the 
dwarf galaxies. The UFD symbols are the same as in Fig. 1211 while 
the dSph data is from this paper and Tafelmeyer et al. (2010). 



of these elements. A neutron-rich, a-freezeout in hyper- 
novae would also contribute to Na production, and thus 
a lack of hypernovae (i.e., if this is the source of the Sr 
excess, as in Farouqi et al. 's 2009 models) could explain 
the very low [Na/Fe] ratios in the most metal-poor stars. 
Evolution models (e.g., Salvadori et al. 2008, Brooks et 
al. 2009, 2007, Ferrara & Tolstoy 2000) suggest that very 
low mass galaxies can lose gas after the first star form- 
ing epoch through SN II driven winds (possibly reac- 
creting cold gas for later star formation events). Koch et 
al. (2006) also suggested that the loss of metal rich winds 
would help to explain the MDF of Carina, which suffers 
from the well known G dwarf problem. Detailed models 
of the star formation and chemical evolution of Carina 
by Lanfranchi & Matteucci (2004) invoked two major 
epochs of outflows through winds. Their best model in- 
cluded a wind efficiency that was 7x the star formation 
rate and associated with the initial and lA star forma- 
tion episodes. Thus, the detailed chemical composition 
of the most metal-poor stars in Carina suggest a lack of 
hypernovae contributions, possibly implying that these 
stars did not form, but more likely indicating that their 
gas was removed from Carina by SN driven winds. 

One alternative to this scenario is that the metal-poor 
stars are displaying the imprint of the hypernovae, rather 
than the lack of hypernovae, which would provide a very 
strong constraint on the nucleosynthetic models. An ex- 
cellent test of this alternative would be the [Zn/Fe] abun- 
dances, which are predicted to be enriched by hypernovae 
because of an increase in the mass ratio between the 
complete and incomplete Si-burning regions (Nomoto et 
al. 2011). We are only able to determine upper limits to 
the [Zn/Fe] ratios in three stars in Carina below [Fe/H] 
= — 2 from our moderate SNR Magellan/MIKE spectra; 
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however these hmits are tantahzingly close to providing 
an interesting constraint (see Fig. [Til) . In particular, Car- 
5070 has a low [Zn/Fe] upper limit, and suggests that 
Car-5070 lacks hypernova enrichments. 

5.3.4. Summary 

Carina and Draco (and possibly Bootes I) appear to be 
quite different in their early chemical evolution from the 
other dSphs and UFDs. These two (three) galaxies may 
be at the critical mass where SN driven winds remove 
the gas from the most massive or energetic SNe II pro- 
genitors, but the products of the remaining SNe II are 
retained, and contribute to the (inhomogeneous) chemi- 
cal evolution of the host. This is unlike the dSphs, which 
appear to retain the gas from the earliest epochs and 
undergo a smooth chemical evolution that is not too dif- 
ferent from stars in the Galactic halo (other than con- 
tributions occurring at lower metallicities, e.g., the AGB 
yields are from metal-poor stars, and the SNe la con- 
tribute at lower metallicities). Both of these galactic sys- 
tems are unlike the UFDs, where a single massive SN II 
may remove all of the gas, quenching the star formation 
event, and imprinting their unique chemical signatures 
on the few stars that will complete the star formation 
process - in this case, any abundance variations within 
an UFD are due to stochastic sampling of that hypernova 
and not chemical evolution. 

5.4. Car-612: A pocket of SN la Enriched Gas 

As discussed in Section 4, Car-612 is underabundant 
in nearly every element when examined relative to iron, 
i.e., [X/Fe]. Therefore, we propose that this star is iron 
enhanced, most likely due to an excess of SN la contri- 
butions in the gas cloud from which it formed. 

In Fig. [26l we show the abundance distribution rela- 
tive to the Galactic averages at [Fe/H] = -1.3 (estimated 
from the Figures in Section 4). Notice that very few ele- 
ment ratios lay in the grey band that would describe the 
[X/Fe] ratios of a typical Galactic star at that metallicity. 
The only exceptions are Cr, Mn, and Fe, and possibly Si 
(though the error is large for that element). If this star 
has an excess of iron by 0.7 dex (a factor of 5), then re- 
moving this in the [X/Fe] ratios would produce the values 
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Fig. 26. — Abundance distribution for Car-612 relative to Galac- 
tic stars at its metallicity of [Fe/H] = —1.3. The Galactic abun- 
dances are estimated from the distributions in Figs. [5] to 1171 The 
grey shaded area represents a (mean) error in the Galactic abun- 
dance of ±0.2 dex. The cyan shaded area represents an offset by 
—0.7 dex, i.e., if Car-612 were to have a metallicity of —2.0. Note 
that at [Fe/H] = —2.0, the stars in classical dSphs have similar 
[X/Fe] values as stars in the Galaxy, thus the apparent iron enrich- 
ment does not depend on differences in the chemical evolution of 
the comparison dataset. 



in the blue band (other than Fe itself which we do not 
adjust). Now Cr, Mn, and Fe (itself) appear enhanced, 
however nearly all other [X/Fe] ratios are in agreement 
with the Galactic stars. A direct comparison of spectral 
line strengths shows that Car-612 has similar Fe I line 
strengths to stars with similar atmospheric parameters 
(e.g., Car-705 in Fig. [Tni also to Car-1013 not shown), 
while the other lines are weaker. 

This peculiar chemistry makes Car-612 similar to three 
iron-rich stars in the outer halo, studied by Ivans et 
al. (2003). These stars are all near [Fe/H] = —2 and have 
low [a/Fe] ratios, low ratios of Y, Sr, and Ba, and two 
show enhancements in Cr, Mn, Ni, and Zn. Of those two 
stars, one is enhanced in Si and Eu. Ivans et al. concluded 
that these stars have larger SN la/II contributions, by 
factors of 3 to 7 relative to the average halo star. This is 
similar to Car-612 where we suggest the enhancement is 
by a factor of 5. 

A deep and thorough examination of the predicted 
yields from existing SN models showed that no combi- 
nation could reproduce the detailed abundance patterns 
of the three outer halo stars, i.e., problems remained in 
the abundances of Ti, Cr, Mn, Ni, and Zn. Since the 
SN yields are able to reproduce the chemistry of the ma- 
jority of stars in the Galaxy, Ivans et al. suggested that 
perhaps the SN yields should not be integrated over a 
Salpeter IMF for these three stars and/or the degree of 
mixing in the ISM may vary from region to region. If 
so, they considered that it is also possible for these stars 
to have been deposited in the outer halo during a dwarf 
galaxy merger. Variations in their detailed abundances 
are therefore related to differences in the chemical evolu- 
tion of their hosts. Car-612 fits with this hypothesis — it 
appears to be iron-enhanced, and yet its detailed chem- 
istry is again different from that of the three outer halo 
stars, e.g., Ni, Zn, and Eu are not enhanced. Identifying 
this star within the Carina dwarf galaxy provides ad- 
ditional evidence for inhomogeneous mixing in this low 
mass galaxy, and provides a clear connection between 
the formation of these chemically peculiar stars in dwarf 
galaxies and their existence in the outer Galactic halo. 

6. CONCLUSIONS 

Carina is an interesting galaxy for chemical evolution 
studies because of its low mass and its episodic star 
formation history. In this paper, we have determined 
the abundances of 23 elements in the spectra of nine 
RGB stars in the Carina dSph galaxy taken with both 
the VLT/FLAMES-UVES and Magellan/MIKE spectro- 
graphs. This is a significant increase from the previous 
number of stars with detailed analyses, e.g., by Shetrone 
et al. (2003, where all stars had intermediate metallic- 
ities) and Koch et al. (2008, where only the iron and 
a-element abundances were determined). Our analysis 
uses both photometric and spectroscopic techniques to 
determine the stellar parameters, and we use new spher- 
ical models, an expanded line list, continuum scattering 
corrections, and hyperfine structure and NLTE correc- 
tions (when available) to improve the precision in our 
abundances. Adopting the ages determined by Lemasle 
et al. (2012), we are able to examine the chemical evolu- 
tion of Carina, separating chemical signatures in the old 
and intermediate- aged populations. A summary of the 
most important results in this paper are as follows: 
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• Inhomogeneous mixing: A largo dispersion in 
[Mg/Fe] indicates poor mixing in the old popula- 
tion. An offset in the [a/Fe] ratios between the 
old and intermediate-aged populations (when pre- 
viously published data are included) also suggests 
that the second star formation event occurred in a- 
enriched gas. In addition, one star, Car-612, seems 
to have formed in a pocket enhanced in SN la/II 
products. 

• SN driven winds: Stars with [Fe/H] < —2 do not 
appear to have been enriched in AGB or SNe la 
products. Their peculiar chemistry includes very 
low ratios of [Sr/Ba], an element ratio that usu- 
ally shows an excess in metal-poor stars in the 
Galaxy and dSphs. Adopting a scenario where 
the excess Sr forms in the more massive or ener- 
getic SN II, then the lack of this excess in Ca- 
rina (also Draco and Bootes I) is consistent with 
the loss of those products by SN II driven winds. 
Low ratios of [Na/Fe], [Mn/Fe], and [Cr/Fe] sup- 
port this scenario, with additional evidence from 
the low [Zn/Fe] upper Hmit for one star, Car-5070. 
The a-elements ratios in Car-5070 are also lower 
than the Galactic distribution. 

It is interesting that the chemistry of the metal-poor 
stars in Carina arc not similar to those in the Galaxy, 
most of the other dSphs, or the UFDs. The [Sr/Fe] and 
[Sr/Ba] ratios are clear indicators of the differences in the 



early chemical evolution of these systems: the Galaxy 
and dSphs appear to retain all of their SN II prod- 
ucts, and show excesses in the [Sr/Ba] ratios, whereas 
Carina (with Draco and Bootes I) may be at the crit- 
ical mass where some gas is lost through SN II driven 
winds, by showing very low [Sr/Fe] and [Sr/Ba] ratios. 
In the UFDs, all of the gas may be lost with the first SN 
II, quenching star formation, and imprinting the unique 
chemical signature of that SN II on the remaining stars 
as they complete the formation process, thus [Sr/Fe] is 
low, but [Sr/Ba] can vary. 

It is also interesting to find a star with an enhancement 
in the SN la/II products in Carina, similar to the three 
outer halo stars examined by Ivans et al. (2003). This 
provides the first direct link between the formation of 
these stars in low mass galaxies and their presence in the 
outer Galactic halo 
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